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Abstract 

In these lectures I review the present status of the so-called Standard Cos- 
mological Model, based on the hot Big Bang Theory and the Inflationary 
Paradigm. I will make special emphasis on the recent developments in obser- 
O ! vational cosmology, mainly the acceleration of the universe, the precise mea- 

surements of the microwave background anisotropics, and the formation of 

■ structure like galaxies and clusters of galaxies from tiny primodial fluctuations 

, generated during inflation. 

' 1. INTRODUCTION 



> 



The last five years have seen the coming of age of Modern Cosmology, a mature branch of science 
based on the hot Big Bang theory and the Inflationary Paradigm. In particular, we can now define rather 
precisely a Standard Model of Cosmology, where the basic parameters are determined within small 
uncertainties, of just a few percent, thanks to a host of experiments and observations. This precision era 
^ \ of cosmology has become possible thanks to important experimental developments in all fronts, from 

O measurements of supernovae at high redshifts to the microwave background anisotropics, as well as to 

^ ■ the distribution of matter in galaxies and clusters of galaxies. 

In these lecture notes I will first introduce the basic concepts and equations associated with hot 
Big Bang cosmology, defining the main cosmological parameters and their corresponding relationships. 
Then I will address in detail the three fundamental observations that have shaped our present knowledge: 
the recent acceleration of the universe, the distribution of matter on large scales and the anisotropics in 
the microwave background. Together these observations allow the precise determination of a handful of 
cosmological parameters, in the context of the inflationary plus cold dark matter paradigm. 

X ' 2. BIG BANG COSMOLOGY 

b : 

Our present understanding of the universe is based upon the successful hot Big Bang theory, which 
explains its evolution from the first fraction of a second to our present age, around 13.6 billion years 
later. This theory rests upon four robust pillars, a theoretical framework based on general relativity, 
as put forward by Albert Einstein Q and Alexander A. Friedmann J2) in the 1920s, and three basic 
observational facts: First, the expansion of the universe, discovered by Edwin P. Hubble in the 1930s, 
as a recession of galaxies at a speed proportional to their distance from us. Second, the relative abundance 
of light elements, explained by George Gamow [4] in the 1940s, mainly that of helium, deuterium and 
lithium, which were cooked from the nuclear reactions that took place at around a second to a few minutes 
after the Big Bang, when the universe was a few times hotter than the core of the sun. Third, the cosmic 
microwave background (CMB), the afterglow of the Big Bang, discovered in 1965 by Arno A. Penzias 
and Robert W. Wilson [ 5 1 as a very isotropic blackbody radiation at a temperature of about 3 degrees 
Kelvin, emitted when the universe was cold enough to form neutral atoms, and photons decoupled from 
matter, approximately 380,000 years after the Big Bang. Today, these observations are confirmed to 
within a few percent accuracy, and have helped establish the hot Big Bang as the preferred model of the 
universe. 

Modern Cosmology begun as a quantitative science with the advent of Einstein's general rela- 
tivity and the realization that the geometry of space-time, and thus the general attraction of matter, is 



determined by the energy content of the universe [ 6 ] 

G^v = R^v - ^g^R + Ag, lu = 8ttG . (1) 

These non-linear equations are simply too difficult to solve without invoking some symmetries of the 
problem at hand: the universe itself. 

We live on Earth, just 8 light-minutes away from our star, the Sun, which is orbiting at 8.5 kpc 
from the center of our galaxy, 1 the Milky Way, an ordinary galaxy within the Virgo cluster, of size a 
few Mpc, itself part of a supercluster of size a few 100 Mpc, within the visible universe, approximately 
10,000 Mpc in size. Although at small scales the universe looks very inhomogeneous and anisotropic, 
the deepest galaxy catalogs like 2dF GRS and SDSS suggest that the universe on large scales (beyond the 
supercluster scales) is very homogeneous and isotropic. Moreover, the cosmic microwave background, 
which contains information about the early universe, indicates that the deviations from homogeneity and 
isotropy were just a few parts per million at the time of photon decoupling. Therefore, we can safely 
impose those symmetries to the univerge at large and determine the corresponding evolution equations. 
The most general metric satisfying homogeneity and isotropy is the Friedmann-Robertson- Walker (FRW) 
metric, written here in terms of the invariant geodesic distance ds 2 = g liu dx^ l dx v in four dimensions [6| 
/i = 0,l,2,3, 2 
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characterized by just two quantities, a scale factor a(t), which determines the physical size of the uni- 
verse, and a constant K, which characterizes the spatial curvature of the universe, 
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Spatially open, flat and closed universes have different three-geometries. Light geodesies on these uni- 
verses behave differently, and thus could in principle be distinguished observationally, as we shall discuss 
later. Apart from the three-dimensional spatial curvature, we can also compute a four-dimensional space- 
time curvature, 

(% = 6°+6^ 2 + 6^. (4) 
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Depending on the dynamics (and thus on the matter/energy content) of the universe, we will have dif- 
ferent possible outcomes of its evolution. The universe may expand for ever, recollapse in the future or 
approach an asymptotic state in between. 



2.1 The matter and energy content of the universe 

The most general matter fluid consistent with the assumption of homogeneity and isotropy is a perfect 
fluid, one in which an observer comoving with the fluid would see the universe around it as isotropic. The 
energy momentum tensor associated with such a fluid can be written as HO 

=pg^ + (p + p) U^U U , (5) 

where p(t) and p(t) are the pressure and energy density of the fluid at a given time in the expansion, as 
measured by this comoving observer, and is the comoving four- velocity, satisfying U^U^ = — 1. For 
such a comoving observer, the matter content looks isotropic (in its rest frame), 

r"„ = d iag(-p(t), p(t), P (t), p(t)) . (6) 

'One parallax second (1 pc), parsec for short, corresponds to a distance of about 3.26 light-years or 3.09 x 10 18 cm. 
2 I am using c = 1 everywhere, unless specified, and a metric signature ( — , +, +, +)• 



The conservation of energy (T fJ-u . u = 0), a direct consequence of the general covariance of the theory 
(G^ v . v = 0), can be written in terms of the FRW metric and the perfect fluid tensor (J5} as 

p + 3-(p + p) = 0. (7) 
a 

In order to find explicit solutions, one has to supplement the conservation equation with an equa- 
tion of state relating the pressure and the density of the fluid, p = p(p). The most relevant fluids in 
cosmology are barotropic, i.e. fluids whose pressure is linearly proportional to the density, p = w p, and 
therefore the speed of sound is constant in those fluids. 

We will restrict ourselves in these lectures to three main types of barotropic fluids: 

• Radiation, with equation of state pn = pn/3, associated with relativistic degrees of freedom (i.e. 
particles with temperatures much greater than their mass). In this case, the energy density of 
radiation decays as pr ~ a -4 with the expansion of the universe. 

• Matter, with equation of state pm — 0, associated with nonrelativistic degrees of freedom (i.e. 
particles with temperatures much smaller than their mass). In this case, the energy density of 
matter decays as pu ~ a~ 3 with the expansion of the universe. 

• Vacuum energy, with equation of state py = —py, associated with quantum vacuum fluctuations. 
In this case, the vacuum energy density remains constant with the expansion of the universe. 
This is all we need in order to solve the Einstein equations. Let us now write the equations of 

motion of observers comoving with such a fluid in an expanding universe. According to general relativity, 
these equations can be deduced from the Einstein equations (Q, by substituting the FRW metric © and 
the perfect fluid tensor ©. The p = i, v = j component of the Einstein equations, together with the 
p = 0, v = component constitute the so-called Friedmann equations, 

a\ 2 8vrG A K 

- = ^P+T"-' (8) 
a J 6 6 o z 
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- = --5-0> + 3p) + T ■ (9) 
a 6 6 

These equations contain all the relevant dynamics, since the energy conservation equation © can be 

obtained from these. 

2.2 The Cosmological Parameters 

I will now define the most important cosmological parameters. Perhaps the best known is the Hubble 
parameter or rate of expansion today, Hq = a/a(to). We can write the Hubble parameter in units of 100 
kms _1 Mpc _1 , which can be used to estimate the order of magnitude for the present size and age of the 
universe, 

H = WOh kms^Mpc -1 , (10) 
cHq 1 = 3000/i _1 Mpc, (11) 
Hq 1 = 9.773 h- 1 Gyr . (12) 

The parameter h was measured to be in the range 0.4 < h < 1 for decades, and only in the last few years 
has it been found to lie within 4% of h = 0.70. I will discuss those recent measurements in the next 
Section. 

Using the present rate of expansion, one can define a critical density p c , that which corresponds 
to a flat universe, 

p c = |^ = I.88/1 2 10~ 29 g/cm 3 (13) 



8vrG 

= 2.77 k- 1 10 11 A/q/^ 1 Mpc) 3 (14) 
= 11.26 h 2 protons/m 3 , (15) 



where M & = 1.989 x 10 33 g is a solar mass unit. The critical density p c corresponds to approximately 
6 protons per cubic meter, certainly a very dilute fluid! 

In terms of the critical density it is possible to define the density parameter 

^|Sp(io) = -(to), (16) 

whose sign can be used to determine the spatial (three-)curvature. Closed universes (K = +1) have 
Qo > 1, flat universes (K = 0) have S7o = 1, and open universes (K = —1) have f^o < 1> no matter 
what are the individual components that sum up to the density parameter. 

In particular, we can define the individual ratios Q i = pi/p c , for matter, radiation, cosmological 
constant and even curvature, today, 

8ttGpm n 8irGp R 

^ A = ^2 n K = -~2TT2- (18) 

For instance, we can evaluate today the radiation component J1r, corresponding to relativistic parti- 
cles, from the density of microwave background photons, pcmb = 7r 2 lc 4 r^ MB /(15?i 3 c 3 ) = 4.5 x 
10 -34 g/cm 3 , which gives ricMB = 2.4 x 10~ 5 h~ 2 . Three approximately massless neutrinos would 
contribute a similar amount. Therefore, we can safely neglect the contribution of relativistic particles to 
the total density of the universe today, which is dominated either by non-relativistic particles (baryons, 
dark matter or massive neutrinos) or by a cosmological constant, and write the rate of expansion in terms 
of its value today, as 

H\a) = Hi (n R 4 + n M 4 + o A + n K a \ ). ( i9) 

\ Ob (X Ob I 

An interesting consequence of these definitions is that one can now write the Friedmann equation today, 
a = ao, as a cosmic sum rule, 

1 = n M + Q A + Q K , (20) 

where we have neglected today. That is, in the context of a FRW universe, the total fraction of 
matter density, cosmological constant and spatial curvature today must add up to one. For instance, if we 
measure one of the three components, say the spatial curvature, we can deduce the sum of the other two. 

Looking now at the second Friedmann equation Q, we can define another basic parameter, the 
deceleration parameter, 

a a AttG r i 

Q0 = —^(to) = [PW + 3 H*o)J , (21) 

defined so that it is positive for ordinary matter and radiation, expressing the fact that the universe expan- 
sion should slow down due to the gravitational attraction of matter. We can write this parameter using 
the definitions of the density parameter for known and unknown fluids (with density Q, x and arbitrary 
equation of state w x ) as 



X 



Uniform expansion corresponds to go = and requires a cancellation between the matter and vacuum 
energies. For matter domination, qo > 0, while for vacuum domination, qo < 0. As we will see in a 
moment, we are at present probing the time dependence of the deceleration parameter and can determine 
with some accuracy the moment at which the universe went from a decelerating phase, dominated by 
dark matter, into an acceleration phase at present, which seems to indicate the dominance of some kind 
of vacuum energy. 




Fig. 1: Parameter space (Q,m, Qa). The green (dashed) line Qa = 1 — Qm corresponds to a flat universe, Qk = 0, separating 
open from closed universes. The blue (dotted) line Qa = corresponds to uniform expansion, go = 0, separating 

accelerating from decelerating universes. The violet (dot-dashed) line corresponds to critical universes, separating eternal 
expansion from recollapse in the future. Finally, the red (continuous) lines correspond to toHo = 0.5, 0.6, . . . , oo, beyond 
which the universe has a bounce. 

2.3 The (fijtf, Q A ) Plane 

Now that we know that the universe is accelerating, one can parametrize the matter/energy content of 
the universe with just two components: the matter, characterized by Q^, and the vacuum energy Qa- 
Different values of these two parameters completely specify the universe evolution. It is thus natural to 
plot the results of observations in the plane (Q m , Qa), hi order to check whether we arrive at a consistent 
picture of the present universe from several different angles (different sets of cosmological observations). 

Moreover, different regions of this plane specify different behaviors of the universe. The bound- 
aries between regions are well defined curves that can be computed for a given model. I will now describe 
the various regions and boundaries. 

• Uniform expansion (qo = 0). Corresponds to the line Qa = Qm/2- Points above this line 
correspond to universes that are accelerating today, while those below correspond to decelerat- 
ing universes, in particular the old cosmological model of Einstein-de Sitter (EdS), with Qa = 
0, Qm = 1. Since 1998, all the data from Supernovae of type la appear above this line, many 
standard deviations away from EdS universes. 

• Flat universe (Qj^ = 0). Corresponds to the line Qa = 1 — Qa/- Points to the right of this line 
correspond to closed universes, while those to the left correspond to open ones. In the last few 
years we have mounting evidence that the universe is spatially flat (in fact Euclidean). 

• Bounce (t^Ho = oo)- Corresponds to a complicated function of Qa(Qm), normally expressed as 
an integral equation, where 

t H = C da [1 + fijif (1/a - 1) + Q A (a 2 - 1)]~ 1/2 
Jo 



is the product of the age of the universe and the present rate of expansion. Points above this line 
correspond to universes that have contracted in the past and have later rebounced. At present, these 
universes are ruled out by observations of galaxies and quasars at high redshift (up to z = 10). 
• Critical Universe (H = H = 0). Corresponds to the boundary between eternal expansion in the 
future and recollapse. For Qm < 1> it is simply the line Q,\ = 0, but for Qm > 1, it is a more 
complicated curve, 

n A = 4n M sm -arcsm — ~— . 

V \Im /j "M 

These critical solutions are asymptotic to the EdS model. 

These boundaries, and the regions they delimit, can be seen in Fig. 1, together with the lines of 
equal toHo values. 

In summary, the basic cosmological parameters that are now been hunted by a host of cosmological 
observations are the following: the present rate of expansion Hq; the age of the universe to; the deceler- 
ation parameter qo; the spatial curvature Q,k\ the matter content QmI the vacuum energy Q^', the baryon 
density f2g; the neutrino density Q u , and many other that characterize the perturbations responsible for 
the large scale structure (LSS) and the CMB anisotropics. 



2.4 The accelerating universe 

Let us first describe the effect that the expansion of the universe has on the objects that live in it. In the 
absence of other forces but those of gravity, the trajectory of a particle is given by general relativity in 
terms of the geodesic equation 

^- + u"u A = , (23) 

where = (7, 7i> 1 ), with 7 2 = 1 — v 2 and v l is the peculiar velocity. Here T^ A is the Christoffel con- 
nection [6], whose only non-zero component is r°- = (a/a) giji substituting into the geodesic equation, 
we obtain \u\ oc 1/a, and thus the particle's momentum decays with the expansion like p oc 1/a. In the 
case of a photon, satisfying the de Broglie relation p = h/X, one obtains the well known photon redshift 

Ai Obit] ) An — Ai On 

-i = -)-i( => z = J^^ = JL_i, (24 ) 
Ao a (*o) Ai ai 

where Ao is the wavelength measured by an observer at time to, while Ai is the wavelength emitted 
when the universe was younger (t\ < to). Normally we measure light from stars in distant galaxies and 
compare their observed spectra with our laboratory (restframe) spectra. The fraction d24l > then gives the 
redshift z of the object. We are assuming, of course, that both the emitted and the restframe spectra are 
identical, so that we can actually measure the effect of the intervening expansion, i.e. the growth of the 
scale factor from t\ to to, when we compare the two spectra. Note that if the emitting galaxy and our own 
participated in the expansion, i.e. if our measuring rods (our rulers) also expanded with the universe, we 
would see no effect! The reason we can measure the redshift of light from a distant galaxy is because our 
galaxy is a gravitationally bounded object that has decoupled from the expansion of the universe. It is 
the distance between galaxies that changes with time, not the sizes of galaxies, nor the local measuring 
rods. 

We can now evaluate the relationship between physical distance and redshift as a function of the 
rate of expansion of the universe. Because of homogeneity we can always choose our position to be at 
the origin r = of our spatial section. Imagine an object (a star) emitting light at time t\, at coordinate 
distance r\ from the origin. Because of isotropy we can ignore the angular coordinates (9, <fi). Then the 
physical distance, to first order, will be d = an n- Since light travels along null geodesies HO, we can 



write = — dt 2 + a 2 (t) dr 2 / (1 — Kr 2 ), and therefore, 



1 f(n) = { n K = o (25) 
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If we now Taylor expand the scale factor to first order, 

1 o(t) 



1 + z ao 
we find, to first approximation, 



1 + H Q (t - t ) + 0(t - to) 2 , (26) 



ri « /(ri) = —(to - ti) + • . . = ~4r + • • • 
a ao-no 

Putting all together we find the famous Hubble law 

Hq d = aoHor\ = z ~ vc, (27) 

which is just a kinematical effect (we have not included yet any dynamics, i.e. the matter content of 
the universe). Note that at low redshift (z <C 1), one is tempted to associate the observed change in 
wavelength with a Doppler effect due to a hypothetical recession velocity of the distant galaxy. This 
is only an approximation. In fact, the redshift cannot be ascribed to the relative velocity of the distant 
galaxy because in general relativity (i.e. in curved spacetimes) one cannot compare velocities through 
parallel transport, since the value depends on the path! If the distance to the galaxy is small, i.e. z <C 1, 
the physical spacetime is not very different from Minkowsky and such a comparison is approximately 
valid. As z becomes of order one, such a relation is manifestly false: galaxies cannot travel at speeds 
greater than the speed of light; it is the stretching of spacetime which is responsible for the observed 
redshift. 



B Band 




Fig. 2: The Type la supernovae observed nearby show a relationship between their absolute luminosity and the timescale of 
their light curve: the brighter supernovae are slower and the fainter ones are faster. A simple linear relation between the absolute 
magnitude and a "stretch factor" multiplying the light curve timescale fits the data quite well. From Ref. (7). 



Hubble's law has been confirmed by observations ever since the 1920s, with increasing precision, 
which have allowed cosmologists to determine the Hubble parameter Ho with less and less systematic 
errors. Nowadays, the best determination of the Hubble parameter was made by the Hubble Space 
Telescope Key Project [SO, Ho = 72 ± 8 km/s/Mpc. This determination is based on objects at distances 
up to 500 Mpc, corresponding to redshifts z < 0.1 . 



Nowadays, we are beginning to probe much greater distances, corresponding to z ~ 1, thanks to 
type la supernovae. These are white dwarf stars at the end of their life cycle that accrete matter from 
a companion until they become unstable and violently explode in a natural thermonuclear explosion 
that out-shines their progenitor galaxy. The intensity of the distant flash varies in time, it takes about 
three weeks to reach its maximum brightness and then it declines over a period of months. Although 
the maximum luminosity varies from one supernova to another, depending on their original mass, their 
environment, etc., there is a pattern: brighter explosions last longer than fainter ones. By studying the 
characteristic light curves, see Fig. 2, of a reasonably large statistical sample, cosmologists from the 
Supernova Cosmology Project Q and the High-redshift Supernova Project [9], are now quite confident 
that they can use this type of supernova as a standard candle. Since the light coming from some of these 
rare explosions has travelled a large fraction of the size of the universe, one expects to be able to infer 
from their distribution the spatial curvature and the rate of expansion of the universe. 



Supernova cosmology project 
Knopetal. (2003) 
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Fig. 3: Upper panel: The Hubble diagram in linear redshift scale. Supernovae with Az < 0.01 of eachother have been 
weighted-averaged binned. The solid curve represents the best-fit flat universe model, (Qm = 0.25, Qa = 0.75). Two other 
cosmological models are shown for comparison, (Qm = 0.25, Qa = 0) and (Qm = 1, Qa = 0). Lower panel: Residuals of 
the averaged data relative to an empty universe. From Ref. 0. 



The connection between observations of high redshift supernovae and cosmological parameters 
is done via the luminosity distance, defined as the distance at which a source of absolute luminosity 
(energy emitted per unit time) C gives a flux (measured energy per unit time and unit area of the detector) 
T = C/Att d\. One can then evaluate, within a given cosmological model, the expression for di as a 



function of redshift [10], 



H d L {z) = iq 11/2 smn 



o v/(l + z') 2 (l + z'Q M ) - z'{2 + z')^A 



(28) 



where sinn(x) = x if K = 0; sin(x) if K = +1 and sinh(a;) if A" = — 1, and we have used the cosmic 
sum rule d20b . 

Astronomers measure the relative luminosity of a distant object in terms of what they call the 
effective magnitude, which has a peculiar relation with distance, 

m(z) = M + 5 log 10 + 25 = M + 5 log 10 [H d L (z)} . (29) 

L Mpc J 

Since 1998, several groups have obtained serious evidence that high redshift supernovae appear fainter 
than expected for either an open (Qm < 1) or a flat (Om = 1) universe, see Fig. 3. In fact, the universe 
appears to be accelerating instead of decelerating, as was expected from the general attraction of matter, 
see Eq. d22l : something seems to be acting as a repulsive force on very large scales. The most natural 
explanation for this is the presence of a cosmological constant, a diffuse vacuum energy that permeates 
all space and, as explained above, gives the universe an acceleration that tends to separate gravitationally 
bound systems from each other. The best-fit results from the Supernova Cosmology Project [ 11 ] give a 
linear combination 

0.8n M - 0.6 U A = -0.16 ±0.05 (ler), 

which is now many sigma away from an EdS model with A = 0. In particular, for a flat universe this 
gives 

n A = 0.71 ±0.05 and Q M = 0.29 ± 0.05 (la). 

Surprising as it may seem, arguments for a significant dark energy component of the universe where 
proposed long before these observations, in order to accommodate the ages of globular clusters, as well 
as a flat universe with a matter content below critical, which was needed in order to explain the observed 
distribution of galaxies, clusters and voids. 

Taylor expanding the scale factor to third order, 

?®=1 + H (t- t ) - ^H*(t - t ) 2 + |# 3 (t - t ) 3 + 0{t - to) 4 , (30) 



where 

qo = ~ ^2 (*o) = g E(! + 3w i) Q i = 2 Um ~ Qa ' (31) 

i 

j = + -4s (*0) = \ E( X + 3W ^( 2 + 3w i) n i = + , (32) 

an 6 2 ~ 

i 

are the deceleration and "jerk" parameters. Substituting into Eq. d28t we find 

#0 d L (z) =z + ^(l-q )z 2 - 1(1 - go - 3q 2 + jo) z 3 + 0(z A ) . (33) 

This expression goes beyond the leading linear term, corresponding to the Hubble law, into the second 
and third order terms, which are sensitive to the cosmological parameters Q, m and Q\. It is only recently 
that cosmological observations have gone far enough back into the early universe that we can begin to 
probe these terms, see Fig. 4. 

This extra component of the critical density would have to resist gravitational collapse, otherwise 
it would have been detected already as part of the energy in the halos of galaxies. However, if most of the 
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Fig. 4: The Supemovae la residual Hubble diagram. Upper panel: Ground-based discoveries are represented by diamonds, 
HST-discovered SNe la are shown as filled circles. Lower panel: The same but with weighted averaged in fixed redshift bins. 
Kinematic models of the expansion history are shown relative to an eternally coasting model q(z) = 0. From Ref. 1 12 . 



energy of the universe resists gravitational collapse, it is impossible for structure in the universe to grow. 
This dilemma can be resolved if the hypothetical dark energy was negligible in the past and only recently 
became the dominant component. According to general relativity, this requires that the dark energy have 
negative pressure, since the ratio of dark energy to matter density goes like a(t)~ 3p ^ p . This argument 
would rule out almost all of the usual suspects, such as cold dark matter, neutrinos, radiation, and kinetic 
energy, since they all have zero or positive pressure. Thus, we expect something like a cosmological 
constant, with a negative pressure, p « —p, to account for the missing energy. 

However, if the universe was dominated by dark matter in the past, in order to form structure, and 
only recently became dominated by dark energy, we must be able to see the effects of the transition from 
the deceleration into the acceleration phase in the luminosity of distant type la supernovae. This has been 
searched for since 1998, when the first convincing results on the present acceleration appeared. However, 
only recently [ 12 1 do we have clear evidence of this transition point in the evolution of the universe. This 
coasting point is defined as the time, or redshift, at which the deceleration parameter vanishes, 

q(z) = -l+(l + z)-^lnH(z)=0, (34) 



where 

H(z) =H \n M (l + zf + n x e^» { ^ Wx(Z )! ^ T +n K (l + zf ~'~ , (35) 



and we have assumed that the dark energy is parametrized by a density Q x today, with a redshift- 
dependent equation of state, w x (z), not necessarily equal to —1. Of course, in the case of a true cosmo- 
logical constant, this reduces to the usual expression. 

Let us suppose for a moment that the barotropic parameter w is constant, then the coasting redshift 
can be determined from 

q[z) 2in M + n x (i + z)to> + n K (i + z)-ii ' K } 



(3\w\ - 1)0, 
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:>,\w\ 



(37) 



which, in the case of a true cosmological constant, reduces to 

Zc = W) - 1 - (38) 

When substituting S1a — 0.7 and Q,m — 0.3, one obtains z c ~ 0.6, in excellent agreement with recent 
observations llT2l . The plane (Qm, ^a) can be seen in Fig. 5, which shows a significant improvement 
with respect to previous data. 




Fig. 5: The recent supernovae data on the (Qm, ^a) plane. Shown are the 1-, 2- and 3-cr contours, as well as the data from 
1998, for comparison. It is clear that the old EdS cosmological model at (SI m = 1, Oa = 0) is many standard deviations away 
from the data. From Ref. [12], 



Now, if we have to live with this vacuum energy, we might as well try to comprehend its origin. 
For the moment it is a complete mystery, perhaps the biggest mystery we have in physics today 1131 . 
We measure its value but we don't understand why it has the value it has. In fact, if we naively predict 
it using the rules of quantum mechanics, we find a number that is many (many!) orders of magnitude 
off the mark. Let us describe this calculation in some detail. In non-gravitational physics, the zero-point 
energy of the system is irrelevant because forces arise from gradients of potential energies. However, we 
know from general relativity that even a constant energy density gravitates. Let us write down the most 
general energy momentum tensor compatible with the symmetries of the metric and that is covariantly 
conserved. This is precisely of the form = Pv 9^ = — Pv 9nv> see Fig- 6. Substituting into 



the Einstein equations Q, we see that the cosmological constant and the vacuum energy are completely 
equivalent, A = 8ttG pv, so we can measure the vacuum energy with the observations of the acceleration 
of the universe, which tells us that — 0.7. 

On the other hand, we can estimate the contribution to the vacuum energy coming from the quan- 
tum mechanical zero-point energy of the quantum oscillators associated with the fluctuations of all quan- 
tum fields, 

$ = E / o A "(0s = ^ E(-D^ + OK 2 4v) , (39) 

where kjjy is the ultraviolet cutoff signaling the scale of new physics. Taking the scale of quantum grav- 
ity, Kjjv = Mpi, as the cutoff, and barring any fortuituous cancellations, then the theoretical expectation 
d39t appears to be 120 orders of magnitude larger than the observed vacuum energy associated with the 
acceleration of the universe, 

pf ~ 1.4 x 10 74 GeV 4 = 3.2 x 10 91 g/cm 3 , (40) 
p^ s ~ 0.7 p c = 0.66 x 1(T 29 g/cm 3 = 2.9 x KT 11 eV 4 . (41) 

Even if we assumed that the ultraviolet cutoff associated with quantum gravity was as low as the elec- 
troweak scale (and thus around the corner, liable to be explored in the LHC), the theoretical expectation 
would still be 60 orders of magnitude too big. This is by far the worst mismatch between theory and 
observations in all of science. There must be something seriously wrong in our present understanding of 
gravity at the most fundamental level. Perhaps we don't understand the vacuum and its energy does not 
gravitate after all, or perhaps we need to impose a new principle (or a symmetry) at the quantum gravity 
level to accommodate such a flagrant mismatch. 



Normal Matter 



a> a> 



d(pV)+pdV = TdS 
A = A 









TdS 


p<0 - 







Vacuum Energy 



Fig. 6: Ordinary matter dilutes as it expands. According to the second law of Thermodynamics, its pressure on the walls should 
be positive, which excerts a force, and energy is lost in the expansion. On the other hand, vacuum energy is always the same, 
independent of the volume of the region, and thus, according to the second law, its pressure must be negative and of the same 
magnitude as the energy density. This negative pressure means that the volume tends to increase more and more rapidly, which 
explains the exponential expansion of the universe dominated by a cosmological constant. 



In the meantime, one can at least parametrize our ignorance by making variations on the idea of a 
constant vacuum energy. Let us assume that it actually evolves slowly with time. In that case, we do not 



expect the equation of state p = —p to remain true, but instead we expect the barotropic parameter w{z) 
to depend on redshift. Such phenomenological models have been proposed, and until recently produced 
results that were compatible with w = — 1 today, but with enough uncertainty to speculate on alternatives 
to a truly constant vacuum energy. However, with the recent supernovae results [ 12 1, there seems to be 
little space for variations, and models of a time-dependent vacuum energy are less and less favoured. In 
the near future, the SNAP satellite [ 14] will measure several thousand supernovae at high redshift and 
therefore map the redshift dependence of both the dark energy density and its equation of state with great 
precision. This will allow a much better determination of the cosmological parameters Qm and 



2.5 Thermodynamics of an expanding plasma 

In this section I will describe the main concepts associated with ensembles of particles in thermal equi- 
librium and the brief periods in which the universe fell out of equilibrium. To begin with, let me make 
contact between the covariant energy conservation law (Q and the second law of thermodynamics, 

TdS = dU + pdV , (42) 

where U = p V is the total energy of the fluid, and p = w p is its barotropic pressure. Taking a comoving 
volume for the universe, V = a 3 , we find 

m dS d , o, d , o. 

T H = H^ + 'li^ = - <43) 

where we have used Q. Therefore, entropy is conserved during the expansion of the universe, dS = 0; 
i.e., the expansion is adiabatic even in those epochs in which the equation of state changes, like in the 
matter-radiation transition (not a proper phase transition). Using Q, we can write 

^ln(pa 3 ) = -3Hw. (44) 

Thus, our universe expands like a gaseous fluid in thermal equilibrium at a temperature T. This tem- 
perature decreases like that of any expanding fluid, in a way that is inversely proportional to the cubic 
root of the volume. This implies that in the past the universe was necessarily denser and hotter. As we 
go back in time we reach higher and higher temperatures, which implies that the mean energy of plasma 
particles is larger and thus certain fundamental reactions are now possible and even common, giving rise 
to processes that today we can only attain in particle physics accelerators. That is the reason why it is 
so important, for the study of early universe, to know the nature of the fundamental interactions at high 
energies, and the basic connection between cosmology and high energy particle physics. However, I 
should clarify a misleading statement that is often used: "high energy particle physics colliders repro- 
duce the early universe" by inducing collisions among relativistic particles. Although the energies of 
some of the interactions at those collisions reach similar values as those attained in the early universe, 
the physical conditions are rather different. The interactions within the detectors of the great particle 
physics accelerators occur typically in the perturbative regime, locally, and very far from thermal equi- 
librium, lasting a minute fraction of a second; on the other hand, the same interactions occurred within 
a hot plasma in equilibrium in the early universe while it was expanding adiabatically and its duration 
could be significantly larger, with a distribution in energy that has nothing to do with those associated 
with particle accelerators. What is true, of course, is that the fundamental parameters corresponding to 
those interactions — masses and couplings — are assumed to be the same, and therefore present terrestrial 
experiments can help us imagine what it could have been like in the early universe, and make predictions 
about the evolution of the universe, in the context of an expanding plasma a high temperatures and high 
densities, and in thermal equilibrium. 



2.51 Fluids in thermal equilibrium 

In order to understand the thermodynamical behaviour of a plasma of different species of particles at high 
temperatures we will consider a gas of particles with g internal degrees of freedom weakly interacting. 
The degrees of freedom corresponding to the different particles can be seen in Table 1. For example, 
leptons and quarks have 4 degrees of freedom since they correspond to the two helicities for both particle 
and antiparticle. However, the nature of neutrinos is still unknown. If they happen to be Majorana 
fermions, then they would be their own antiparticle and the number of degrees of freedom would reduce 
to 2. For photons and gravitons (without mass) their 2 d.o.f. correspond to their states of polarization. 
The 8 gluons (also without mass) are the gauge bosons responsible for the strong interaction betwen 
quarks, and also have 2 d.o.f. each. The vector bosons W ± and Z° are massive and thus, apart from the 
transverse components of the polarization, they also have longitudinal components. 



Particle 


Spin 


Degrees of freedom (g) 


Nature 


Higgs 
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Massive scalar 


photon 


1 


2 


Massless vector 


graviton 


2 


2 


Massless tensor 


gluon 


1 


2 


Massless vector 


Wy Z 


1 


3 


Massive vector 


leptons & quarks 


1/2 


4 


Dirac Fermion 


neutrinos 


1/2 


4(2) 


Dirac (Majorana) Fermion 



Table 1 : The internal degrees of freedom of various fundamental particles. 



For each of these particles we can compute the number density n, the energy density p and the 
pressure p, in thermal equilibrium at a given temperature T, 

n = 5 /(2^ /(p) ' (45) 

p=g /(S^ (p)/(p) ' (46) 



(2tt) 2 

^/(M /(P) ' ^ 



d 3 p |p| 2 

(2vr) 3 3E 

where the energy is given by E 2 = |p| 2 + m? and the momentum distribution in thermal (kinetic) 
equilibrium is 

1 f — 1 Bose — Einstein 

* ( p ) _ e (E-»)/T ± i j +1 Fermi -Dirac ^ 

The chemical potential p is conserved in these reactions if they are in thermal equilibrium. For example, 
for reactions of the type i + j < — ► k + I , we have pi + pj = p^ + p\. For example, the chemical 
potencial of the photon vanishes p 1 = 0, and thus particles and antiparticles have opposite chemical 
potentials. 

From the equilibrium distributions one can obtain the number density n, the energy p and the 
pressure p, of a particle of mass m with chemical potential p at the temperature T, 



P ~2^L e (E- tl )/T± 1 



L 

v= hl 



oo (#2 _ m 2 )3 /2 
m 



e (E-n)/T ± l 

For a non-degenerate <C T) relativistic gas (m <C T), we find 

f C(3) nT 3 

E 2 dE —9? 



(50) 
(51) 
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(52) 



P = - 3 P, 



(53) 



(54) 



where £(3) = 1.20206 ... is the Riemann Zeta function. For relativistic fluids, the energy density per 
particle is 

Bosons 



7T 



11 



30C(3) 

77T 2 



T ~ 2.701 T 
T ~ 3.151 T 



(55) 
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For relativistic bosons or fermions with // < and |/i| < T, we have 



Fermions 
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7T 



1 

^=3 P - 



(56) 
(57) 
(58) 



For a bosonic particle, a positive chemical potential, fi > 0, indicates the presence of a Bose-Einstein 
condensate, and should be treated separately from the rest of the modes. 

On the other hand, for a non-relativistic gas (m » T), with arbitrary chemical potential fj,, we find 



n = g 



p = mn , 
p = nT <C p. 
The average energy density per particle is 



n 



(59) 

(60) 
(61) 

(62) 



Note that, at any given temperature T, the contribution to the energy density of the universe coming 
from non-relativistic particles in thermal equilibrium is exponentially supressed with respect to that of 



relativistic particles, therefore we can write 
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7T 
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g*(T) 



bosons fermions 



(63) 
(64) 



where the factor 7/8 takes into account the difference between the Fermi and Bose statistics; is the 
total number of light d.o.f. (m <C T), and we have also considered the possibility that particle species i 
(bosons or fermions) have an equilibrium distribution at a temperature Ti different from that of photons, 
as happens for example when a given relativistic species decouples from the thermal bath, as we will 
discuss later. This number, , strongly depends on the temperature of the universe, since as it expands 
and cools, different particles go out of equilibrium or become non-relativistic (m S> T) and thus become 
exponentially suppressed from that moment on. A plot of the time evolution of g* (T) can be seen in 
Fig. 7. 
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Fig. 7: the light degrees of freedom g f and g*s as a function of the temperature of the universe. FromRef. [B). 



For example, for T <C 1 MeV, i.e. after the time of primordial Big Bang Nucleosynthesis (BBN) 
and neutrino decoupling, the only relativistic species are the 3 light neutrinos and the photons; since the 
temperature of the neutrinos is T v = (4/H) 1 / 3 ^ = 1.90 K, see below, we have g* = 2 + 3x|x 

(lj) 4/3 = 3 - 36 > while 9*s = 2 + 3 x | x (A) = 3.91. 

For 1 MeV <C T <C 100 MeV, i.e. between BBN and the phase transition from a quark-gluon 
plasma to hadrons and mesons, we have, as relativistic species, apart from neutrinos and photons, also 
the electrons and positrons, sog* = 2 + 3x| + 2x| = 10.75. 

For T S> 250 GeV, i.e. above the electroweak (EW) symmetry breaking scale, we have one photon 
(2 polarizations), 8 gluons (massless), the and Z° (massive), 3 families of quarks & leptones, a Higgs 
(still undiscovered), with which one finds g* = ^ = 106.75. 

At temperatures well above the electroweak transition we ignore the number of d.o.f. of particles, 
since we have never explored those energies in particle physics accelerators. Perhaps in the near future, 
with the results of the Large Hadron Collider (LHC) at CERN, we may may predict the behaviour of the 
universe at those energy scales. For the moment we even ignore whether the universe was in thermal 
equilibrium at those temperatures. The highest energy scale at which we can safely say the universe 



was in thermal equilibrium is that of BBN, i.e. 1 MeV, due to the fact that we observe the present 
relative abundances of the light element produced at that time. For instance, we can't even claim that the 
universe went through the quark-gluon phase transtion, at ~ 200 MeV, since we have not observed yet 
any signature of such an event, not to mention the electroweak phase transition, at ~ 1 TeV. 

Let us now use the relation between the rate of expansion and the temperature of relativistic parti- 
cles to obtain the time scale of the universe as a function of its temperature, 

1/2 T 2 1 _i/ 2 M P f T \- 2 

thus, e.g. at the EW scale (100 GeV) the universe was just 10~ 10 s old, while during the primordial BBN 
(1-0.1 MeV), it was 1 s to 3 min old. 



2.52 The entropy of the universe 

During most of the history of the universe, the rates of reaction, r int , of particles in the thermal bath 
are much bigger than the rate of expansion of the universe, H, so that local thermal equilibrium was 
mantained. In this case, the entropy per comoving volume remained constant. In an expanding universe, 
the second law of thermodynamics, applied to the element of comoving volume, of unit coordinate 
volume and physical volume V = a 3 , can be written as, see d?2l . 

TdS = d(p V)+pdV = d[(p + p)V] - V dp . (66) 

d 2 S d 2 S 

Using the Maxwell condition of integrability, „ mOTr = — — — — , we find that dp = (p + p)dT IT, so that 

oToV oVoT 



dS = d 



(p + p)— + const. 



(67) 



i.e. the entropy in a comoving volume is S = (p + p)^ , except for a constant. Using now the first law, 
the covariant conservation of energy, T^ v v = 0, we have 



d 



(p + p)V]=Vdp =}► d ((p + p)^)=0, (68) 



and thus, in thermal equilibrium, the total entropy in a comoving volume, S = a 3 (p + p)/T, is con- 
served. During most of the evolution of the universe, this entropy was dominated by the contribution 
from relativistic particles, 

2tt 2 

S = — — g* s (aT) 3 = const. , (69) 
45 

9*(T)= £ 9l (^) 3 + l E 9i(^)\ (70) 

bosons fermions 

where g*s is the number of "entropic" degrees of freedom, as we can see in Fig. 7. Above the electron- 
positron annihilation, all relativistic particles had the same temperature and thus g*s = 9*- It may be also 
useful to realize that the entropy density, s = S/a 3 ,is propotional to the number density of relativistic 
particles, and in particular to the number density of photons, s = 1.80g*sn 7 ; today, s = 7.04 n 7 . 
However, since g*s in general is a function of temperature, we can't always interchange s and n 7 . 

The conservation of S implies that the entropy density satisfies s oc a~ 3 , and thus the physical size 
of the comoving volume is a 3 oc s _1 ; therefore, the number of particles of a given species in a comoving 



volume, N = a 3 n, is proportional to the number density of that species over the entropy density s, 



f 45C(3)g 

— T > m, [x 

n~- = \ 27T A g ; s , (7i) 

If this number does not change, i.e. if those particles are neither created nor destroyed, then n/s remains 
constant. As a useful example, we will consider the barionic number in a comoving volume, 

s s 

As long as the interactions that violate barion number occur sufficiently slowly, the barionic number per 
comoving volume, ng/s, will remain constant. Although 

nig Tin 
V= — = l.80g*s — , (73) 

ra 7 s 

the ratio between barion and photon numbers it does not remain constant during the whole evolution 
of the universe since g*$ varies; e.g. during the annihilation of electrons and positrons, the number of 
photons per comoving volume, N y = a 3 n 7 , grows a factor 11/4, and r\ decreases by the same factor. 
After this epoch, however, g* is constant so that rj ~ 7ne/s and tib/s can be used indistinctly. 

Another consequence of Eq, d69l is that 5 = const, implies that the temperature of the universe 
evolves as 

Tong-^a- 1 . (74) 

As long as g*s remains constant, we recover the well known result that the universe cools as it expands 
according to T oc l/a. The factor g^^ appears because when a species becomes non-relativistic (when 
T < m), and effectively disappears from the energy density of the universe, its entropy is transferred 
to the rest of the relativistic particles in the plasma, making T decrease not as quickly as l/a, until g*s 
again becomes constant. 

From the observational fact that the universe expands today one can deduce that in the past it must 
have been hotter and denser, and that in the future it will be colder and more dilute. Since the ratio of 
scale factors is determined by the redshift parameter z, we can obtain (to very good approximation) the 
temperature of the universe in the past with 

T = T (1 + z) . (75) 

This expression has been spectacularly confirmed thanks to the absorption spectra of distant quasars [ 16 1. 
These spectra suggest that the radiation background was acting as a thermal bath for the molecules in the 
interstellar medium with a temperature of 9 K at a redshift z ~ 2, and thus that in the past the photon 
background was hotter than today. Furthermore, observations of the anisotropics in the microwave back- 
ground confirm that the universe at a redshift z = 1089 had a temperature of 0.3 eV, in agreement with 
Eq. G3- 



2.6 The thermal evolution of the universe 

In a strict mathematical sense, it is impossible for the universe to have been always in thermal equilibrium 
since the FRW model does not have a timelike Killing vector. In practice, however, we can say that 
the universe has been most of its history very close to thermal equilibrium. Of course, those periods 
in which there were deviations from thermal equilibrium have been crucial for its evolution thereafter 
(e.g. baryogenesis, QCD transition, primordial nucleosynthesis, recombination, etc.); without these the 
universe today would be very different and probably we would not be here to tell the story. 



The key to understand the thermal history of the universe is the comparison between the rates of 
interaction between particles (microphysics) and the rate of expansion of the universe (macrophysics). 
Ignoring for the moment the dependence of g* on temperature, the rate of change of T is given directly 
by the rate of expansion, T/T = —H. As long as the local interactions — necessary in order that the 
particle distribution function adjusts adiabatically to the change of temperature — are sufficiently fast 
compared with the rate of expansion of the universe, the latter will evolve as a succession of states very 
close to thermal equilibrium, with a temperature proportional to a -1 . If we evaluate the interaction rates 
as 

T int = (na\v\) , (76) 

where n(t) is the number density of target particles, a is the cross section on the interaction and v is the 
relative velocity of the reaction, all averaged on a thermal distribution; then a rule of thumb for ensuring 
that thermal equilibrium is maintained is 

T hlt >H. (77) 

This criterium is understandable. Suppose, as often occurs, that the interaction rate in thermal equilib- 
rium is Tint oc T n , with n > 2; then, the number of interactions of a particle after time t is 

roo i p. 

N int = T int (t')dt' = jj-^t), (78) 

Jt n — 2H 

therefore the particle interacts less than once from the moment in which T mt rs H. If Ti nt ^ H, the 
species remains coupled to the thermal plasma. This doesn't mean that, necessarily, the particle is out 
of local thermal equilibrium, since we have seen already that relativistic particles that have decoupled 
retain their equilibrium distribution, only at a different temperature from that of the rest of the plasma. 

In order to obtain an approximate description of the decoupling of a particle species in an expand- 
ing universe, let us consider two types of interaction: 

i) interactions mediated by massless gauge bosons, like for example the photon. In this case, the cross 
section for particles with significant momentum transfer can be written as a ~ a 2 /T 2 , with a = g 2 /47r 
the coupling constant of the interaction. Assuming local thermal equilibrium, n(t) ~ T 3 and thus the 
interaction rate becomes V ~ n a \ v\ ~ a 2 T. Therefore, 

r 2 M P 

— ~ a — , (79) 

so that for temperatures of the universe T ^ a 2 Mp ~ 10 16 GeV, the reactions are fast enough and the 
plasma is in equilibrium, while for T ^ 10 16 GeV, reactions are too slow to maintain equilibrium and 
it is said that they are "frozen-out". An important consequence of this result is that the universe could 
never have been in thermal equilibrium above the grand unification (GUT) scale. 

ii) interactions mediated by massive gauge bosons, e.g. like the and Z°, or those responsible for the 
GUT interactions, X and Y. We will generically call them X bosons. The cross section depends rather 
strongly on the temperature of the plasma, 

G\T 2 T<^M X 
a ~ { 2 (80) 
^ T » M X 

where Gx ~ a/M x is the effective coupling constant of the interaction at energies well below the mass 
of the vector boson, analogous to the Fermi constant of the electroweak interaction, Gf = g 2 / (Ay/2M^) 
at tree level. Note that for T » Mx we recover the result for massless bosons, so we will concentrate 
here on the other case. For T < Mx, the rate of thermal interactions is T ~ n a \ v\ ~ G 2 X T 5 . Therefore, 

J-~G 2 x MpT 3 , (81) 



such that at temperatures in the range 

M X >T> G x 2/3 M- 1/3 ~ (j^gy) MeV , (82) 

reactions occur so fast that the plasma is in thermal equilibrium, while for T < (M x /100 GeV) 4 / 3 MeV, 
those reactions are too slow for maintaining equilibrium and they effective freeze-out, see Eq. d78b . 



2.61 The decoupling of relativistic particles 

Those relativistic particles that have decoupled from the thermal bath do not participate in the transfer of 
entropy when the temperature of the universe falls below the mass thershold of a given species T ~ m; in 
fact, the temperature of the decoupled relativistic species falls as T oc 1/a, as we will now show. Suppose 
that a relativistic particle is initially in local thermal equilibrium, and that it decoples at a temperature 
Td and time tp- The phase space distribution at the time of decoupling is given by the equilibrium 
distribution, 

f(p,t D ) = 7WJ ^ T - i . (83) 

After decoupling, the energy of each massless particle suffers redshift, E(t) = Erj (au/a(i)). The 
number density of particles also decreases, n(t) = n_o (ai)/a(i)) 3 . Thus, the phase space distribution at 
a time t > t d is 

dPn cl 11 

/(P ' t) = d^= /(P ^' tD) = eW*nT D ± J = ^EjTTi ' (84) 

so that we conclude that the distribution function of a particle that has decoupled while being relativistic 
remains self-similar as the universe expands, with a temperature that decreases as 

T = T D — oca' 1 , (85) 
a 

and not as g~g^ a~ l , like the rest of the plasma in equilibrium (T74b. 



2.62 The decoupling of non-relativistic particles 

Those particles that decoupled from the thermal bath when they were non-relativistic (m 3> T) behave 
differently. Let us study the evolution of the distribution function of a non-relativistic particle that was 
in local thermal equilibrium at a time try, when the universe had a temperature Tjj. The moment of 
each particle suffers redshift as the universe expands, |p| = \pd\ {id/o), see Eq. d24b . Therefore, their 
kinetic energy satisfies E = Ed {cld/cl) 2 . On the other hand, the particle number density also varies, 
n(t) =njj (au/a(t)) 3 , so that a decoupled non-relativistic particle will have an equilibrium distribution 
function characterized by a temperature 

a 2 

T = To — oc a~ 2 , (86) 

and a chemical potential 

T 

H(t) = m + (hd - m) — , (87) 
-Id 

whose variation is precisely that which is needed for the number density of particle to decrease as a~ 3 . 

In summary, a particle species that decouples from the thermal bath follows an equilibrium dis- 
tribution function with a temperature that decreases like Tr oc a -1 for relativistic particles (To 3> m) 
or like Tnr oc a~ 2 for non-relativistic particles (To <C m). On the other hand, for semi-relativistic 
particles (To ~ m), its phase space distribution does not maintain an equilibrium distribution function, 
and should be computed case by case. 



2.63 Brief thermal history of the universe 

I will briefly summarize here the thermal history of the universe, from the Planck era to the present. As 
we go back in time, the universe becomes hotter and hotter and thus the amount of energy available for 
particle interactions increases. As a consequence, the nature of interactions goes from those described at 
low energy by long range gravitational and electromagnetic physics, to atomic physics, nuclear physics, 
all the way to high energy physics at the electroweak scale, gran unification (perhaps), and finally quan- 
tum gravity. The last two are still uncertain since we do not have any experimental evidence for those 
ultra high energy phenomena, and perhaps Nature has followed a different path. 

The way we know about the high energy interactions of matter is via particle accelerators, which 
are unravelling the details of those fundamental interactions as we increase in energy. However, one 
should bear in mind that the physical conditions that take place in our high energy colliders are very 
different from those that occurred in the early universe. These machines could never reproduce the 
conditions of density and pressure in the rapidly expanding thermal plasma of the early universe. Nev- 
ertheless, those experiments are crucial in understanding the nature and rate of the local fundamental 
interactions available at those energies. What interests cosmologists is the statistical and thermal proper- 
ties that such a plasma should have, and the role that causal horizons play in the final outcome of the early 
universe expansion. For instance, of crucial importance is the time at which certain particles decoupled 
from the plasma, i.e. when their interactions were not quick enough compared with the expansion of the 
universe, and they were left out of equilibrium with the plasma. 

One can trace the evolution of the universe from its origin till today. There is still some speculation 
about the physics that took place in the universe above the energy scales probed by present colliders. 
Nevertheless, the overall layout presented here is a plausible and hopefully testable proposal. According 
to the best accepted view, the universe must have originated at the Planck era (10 19 GeV, 10~ 43 s) 
from a quantum gravity fluctuation. Needless to say, we don't have any experimental evidence for such 
a statement: Quantum gravity phenomena are still in the realm of physical speculation. However, it 
is plausible that a primordial era of cosmological inflation originated then. Its consequences will be 
discussed below. Soon after, the universe may have reached the Grand Unified Theories (GUT) era (10 16 
GeV, 10~ 35 s). Quantum fluctuations of the inflaton field most probably left their imprint then as tiny 
perturbations in an otherwise very homogenous patch of the universe. At the end of inflation, the huge 
energy density of the inflaton field was converted into particles, which soon thermalized and became the 
origin of the hot Big Bang as we know it. Such a process is called reheating of the universe. Since 
then, the universe became radiation dominated. It is probable (although by no means certain) that the 
asymmetry between matter and antimatter originated at the same time as the rest of the energy of the 
universe, from the decay of the inflaton. This process is known under the name of baryogenesis since 
baryons (mostly quarks at that time) must have originated then, from the leftovers of their annihilation 
with antibaryons. It is a matter of speculation whether baryogenesis could have occurred at energies 
as low as the electroweak scale (100 GeV, 10~ 10 s). Note that although particle physics experiments 
have reached energies as high as 100 GeV, we still do not have observational evidence that the universe 
actually went through the EW phase transition. If confirmed, baryogenesis would constitute another 
"window" into the early universe. As the universe cooled down, it may have gone through the quark- 
gluon phase transition (10 2 MeV, 10 -5 s), when baryons (mainly protons and neutrons) formed from 
their constituent quarks. 

The furthest window we have on the early universe at the moment is that of primordial nucleosyn- 
thesis (1 — 0.1 MeV, 1 s - 3 min), when protons and neutrons were cold enough that bound systems 
could form, giving rise to the lightest elements, soon after neutrino decoupling: It is the realm of nuclear 
physics. The observed relative abundances of light elements are in agreement with the predictions of the 
hot Big Bang theory. Immediately afterwards, electron-positron annihilation occurs (0.5 MeV, 1 min) 
and all their energy goes into photons. Much later, at about (1 eV, ~ 10 5 yr), matter and radiation have 
equal energy densities. Soon after, electrons become bound to nuclei to form atoms (0.3 eV, 3 x 10 5 



yr), in a process known as recombination: It is the realm of atomic physics. Immediately after, photons 
decouple from the plasma, travelling freely since then. Those are the photons we observe as the cosmic 
microwave background. Much later (~ 1 — 10 Gyr), the small inhomogeneities generated during inflation 
have grown, via gravitational collapse, to become galaxies, clusters of galaxies, and superclusters, char- 
acterizing the epoch of structure formation. It is the realm of long range gravitational physics, perhaps 
dominated by a vacuum energy in the form of a cosmological constant. Finally (3K, 13 Gyr), the Sun, 
the Earth, and biological life originated from previous generations of stars, and from a primordial soup 
of organic compounds, respectively. 

I will now review some of the more robust features of the Hot Big Bang theory of which we have 
precise observational evidence. 

2.64 Primordial nucleosynthesis and light element abundance 

In this subsection I will briefly review Big Bang nucleosynthesis and give the present observational 
constraints on the amount of baryons in the universe. In 1920 Eddington suggested that the sun might 
derive its energy from the fusion of hydrogen into helium. The detailed reactions by which stars burn 
hydrogen were first laid out by Hans Bethe in 1939. Soon afterwards, in 1946, George Gamow realized 
that similar processes might have occurred also in the hot and dense early universe and gave rise to the 
first light elements [4]. These processes could take place when the universe had a temperature of around 
T NS ~ 1 — 0.1 MeV, which is about 100 times the temperature in the core of the Sun, while the density 
is p NS = fo^-^Ns ~ ^ gem -3 , about the same density as the core of the Sun. Note, however, that 
although both processes are driven by identical thermonuclear reactions, the physical conditions in star 
and Big Bang nucleosynthesis are very different. In the former, gravitational collapse heats up the core of 
the star and reactions last for billions of years (except in supernova explosions, which last a few minutes 
and creates all the heavier elements beyond iron), while in the latter the universe expansion cools the hot 
and dense plasma in just a few minutes. Nevertheless, Gamow reasoned that, although the early period of 
cosmic expansion was much shorter than the lifetime of a star, there was a large number of free neutrons 
at that time, so that the lighter elements could be built up quickly by succesive neutron captures, starting 
with the reaction n + p — > D + 7. The abundances of the light elements would then be correlated with 
their neutron capture cross sections, in rough agreement with observations ll6l fT71 . 

Nowadays, Big Bang nucleosynthesis (BBN) codes compute a chain of around 30 coupled nuclear 
reactions [ 18], to produce all the light elements up to beryllium-7. 3 Only the first four or five elements 
can be computed with accuracy better than 1% and compared with cosmological observations. These 
light elements are H, 4 He, D, 3 He, 7 Li, and perhaps also 6 Lz. Their observed relative abundance to 
hydrogen is [1 : 0.25 : 3 • 10~ 5 : 2 • 10~ 5 : 2 • 10 -10 ] with various errors, mainly systematic. The BBN 
codes calculate these abundances using the laboratory measured nuclear reaction rates, the decay rate of 
the neutron, the number of light neutrinos and the homogeneous FRW expansion of the universe, as a 
function of only one variable, the number density fraction of baryons to photons, 77 = ne/w 7 . In fact, 
the present observations are only consistent, see Fig. 8 and Ref. fP71[T8l fl91. with a very narrow range of 
values of 

7/10 = 10 10 77 = 6.2 ± 0.6 . (88) 

Such a small value of 77 indicates that there is about one baryon per 10 9 photons in the universe today. 
Any acceptable theory of baryogenesis should account for such a small number. Furthermore, the present 
baryon fraction of the critical density can be calculated from 7710 as 

n B h 2 = 3.6271 x 10~ 3 7/10 = 0.0224 ± 0.0022 (95% c.l.) (89) 

Clearly, this number is well below closure density, so baryons cannot account for all the matter in the 
universe, as I shall discuss below. 

3 The rest of nuclei, up to iron (Fe), are produced in heavy stars, and beyond Fe in novae and supernovae explosions. 
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Fig. 8: The relative abundance of light elements to Hidrogen. Note the large range of scales involved. From Ref. 1171 . 



2.65 Neutrino decoupling 

Just before the nucleosynthesis of the lightest elements in the early universe, weak interactions were too 
slow to keep neutrinos in thermal equilibrium with the plasma, so they decoupled. We can estimate the 
temperature at which decoupling occurred from the weak interaction cross section, a w ~ G 2 F T 2 at finite 
temperature T, where Gf = 1.2 x 1CP 5 GeV 2 is the Fermi constant. The neutrino interaction rate, via 
W boson exchange in n + v «-> p + e~ and p + v n + e + , can be written as llT31 

T„ = n u (a w \v\)-G F T 5 , (90) 

while the rate of expansion of the universe at that time (g* = 10.75) was H ~ 5.4 T 2 /Mp, where 
Mp = 1.22 x 10 19 GeV is the Planck mass. Neutrinos decouple when their interaction rate is slower 
than the universe expansion, T v < H or, equivalently, at T u _& cc ~ 0.8 MeV. Below this temperature, 
neutrinos are no longer in thermal equilibrium with the rest of the plasma, and their temperature continues 
to decay inversely proportional to the scale factor of the universe. Since neutrinos decoupled before 
e + e~ annihilation, the cosmic background of neutrinos has a temperature today lower than that of the 
microwave background of photons. Let us compute the difference. At temperatures above the the mass 
of the electron, T > m e = 0.511 MeV, and below 0.8 MeV, the only particle species contributing to 
the entropy of the universe are the photons (g* = 2) and the electron-positron pairs (g* = 4 x |); total 
number of degrees of freedom g # = At temperatures T ~ m e , electrons and positrons annihilate into 
photons, heating up the plasma (but not the neutrinos, which had decoupled already). At temperatures 
T < m e , only photons contribute to the entropy of the universe, with g* = 2 degrees of freedom. 



Therefore, from the conservation of entropy, we find that the ratio of T 7 and T u today must be 

T^ = (t) V3 = L401 ^ r„ = 1.945 K, (91) 

where I have used T CMB = 2.725 ± 0.002 K. We still have not measured such a relic background of 
neutrinos, and probably will remain undetected for a long time, since they have an average energy of 
order 10~ 4 eV, much below that required for detection by present experiments (of order GeV), precisely 
because of the relative weakness of the weak interactions. Nevertheless, it would be fascinating if, in the 
future, ingenious experiments were devised to detect such a background, since it would confirm one of 
the most robust features of Big Bang cosmology. 



2.66 Matter-radiation equality 

Relativistic species have energy densities proportional to the quartic power of temperature and therefore 
scale as pr oc a -4 , while non-relativistic particles have essentially zero pressure and scale as pyi oc a~ 3 . 
Therefore, there will be a time in the evolution of the universe in which both energy densities are equal 
PR(*eq) = PM.(teq)- Since then both decay differently, and thus 



1 + Zcq = = = 3.1 x 10 4 n M h 2 , (92) 



where I have used Or/i 2 = Q CMB h 2 + Q, u h 2 = 3.24 x 10 5 for three massless neutrinos at T = T v . As 
I will show later, the matter content of the universe today is below critical, Qm — 0.3, while h ~ 0.71, 
and therefore (1 + z cq ) ~ 3400, or about i eq = 1308 (fi M fr 2 )~V ~ 61,000 years after the origin of 
the universe. Around the time of matter-radiation equality, the rate of expansion (fl9t can be written as 

(o = 1) 

H(a) = H (n R a" 4 + ft M a" 3 ) ^ = H < 2 a~ 3 / 2 (l + ^) V ' . (93) 

The horizon size is the coordinate distance travelled by a photon since the beginning of the universe, 
(in ~ H^ 1 , i.e. the size of causally connected regions in the universe. The comoving horizon size is 
then given by 

Thus the horizon size at matter-radiation equality (a = a eq ) is 

d H (a eq ) = n M 1/2 a^ 2 ~ 12 (ttMh)- 1 h^Mpc . (95) 

This scale plays a very important role in theories of structure formation. 



2.67 Recombination and photon decoupling 

As the temperature of the universe decreased, electrons could eventually become bound to protons to 
form neutral hydrogen. Nevertheless, there is always a non-zero probability that a rare energetic photon 
ionizes hydrogen and produces a free electron. The ionization fraction of electrons in equilibrium with 
the plasma at a given temperature is given by the Saha equation [ 15 ] 

1 - 4V2C(3) / T \ 3 / 2 p /T 

X e H J-k \m P J 



where E- lon = 13.6 eV is the ionization energy of hydrogen, and r\ is the baryon-to-photon ratio dSSl . If 
we now use Eq. ( 1731 . we can compute the ionization fraction X^ q as a function of redshift z. Note that 



the huge number of photons with respect to electrons (in the ratio 4 He : H : 7 ~ 1 : 4 : 10 10 ) implies 
that even at a very low temperature, the photon distribution will contain a sufficiently large number of 
high-energy photons to ionize a significant fraction of hydrogen. In fact, defining recombination as the 
time at which X° q = 0.1, one finds that the recombination temperature is T rec = 0.31 eV <C E- lon , 
for 7710 — 6.2. Comparing with the present temperature of the microwave background, we deduce the 
corresponding redshift at recombination, (1 + z rcc ) ~ 1331. 

Photons remain in thermal equilibrium with the plasma of baryons and electrons through elastic 
Thomson scattering, with cross section 



8vra 2 
3m, 2 



6.65 x 10~ 25 cm 2 



0.665 barn . 



(97) 



where a = 1/137.036 is the dimensionless electromagnetic coupling constant. The mean free path of 
photons A 7 in such a plasma can be estimated from the photon interaction rate, A^ 1 ~ T 7 = n e a T . For 
temperatures above a few eV, the mean free path is much smaller that the causal horizon at that time and 
photons suffer multiple scattering: the plasma is like a dense fog. Photons will decouple from the plasma 
when their interaction rate cannot keep up with the expansion of the universe and the mean free path 
becomes larger than the horizon size: the universe becomes transparent. We can estimate this moment 
by evaluating T 7 = H at photon decoupling. Using n e = X e r]n^, one can compute the decoupling 
temperature as T^ ec = 0.26 eV, and the corresponding redshift as 1 + z^ ec ~ 1100. Recently, WMAP 
measured this redshift to be 1 + z^ cc ~ 1089 ± 1 [20]. This redshift defines the so called last scattering 
surface, when photons last scattered off protons and electrons and travelled freely ever since. This 
decoupling occurred when the universe was approximately t dec = 1.5 x 10 5 (tt M h 2 y 1/2 ~ 380,000 
years old. 
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Fig. 9: The Cosmic Microwave Background Spectrum seen by the FIRAS instrument on COBE. The left panel corresponds to 
the monopole spectrum, To — 2.725 ± 0.002 K, where the error bars are smaller than the line width. The right panel shows the 
dipole spectrum, <5Ti = 3.372 ± 0.014 mK. From Ref. (2TJ. 



2.68 The microwave background 

One of the most remarkable observations ever made my mankind is the detection of the relic background 
of photons from the Big Bang. This background was predicted by George Gamow and collaborators 
in the 1940s, based on the consistency of primordial nucleosynthesis with the observed helium abun- 
dance. They estimated a value of about 10 K, although a somewhat more detailed analysis by Alpher and 
Herman in 1950 predicted T 7 5 K. Unfortunately, they had doubts whether the radiation would have 
survived until the present, and this remarkable prediction slipped into obscurity, until Dicke, Peebles, 



Roll and Wilkinson [22 1 studied the problem again in 1965. Before they could measure the photon back- 
ground, they learned that Penzias and Wilson had observed a weak isotropic background signal at a radio 
wavelength of 7.35 cm, corresponding to a blackbody temperature of T 7 = 3.5 ± 1 K. They published 
their two papers back to back, with that of Dicke et al. explaining the fundamental significance of their 
measurement [6|. 

Since then many different experiments have confirmed the existence of the microwave background. 
The most outstanding one has been the Cosmic Background Explorer (COBE) satellite, whose FIRAS 
instrument measured the photon background with great accuracy over a wide range of frequencies (y = 
1 — 97 cm -1 ), see Ref. [21 ], with a spectral resolution — = 0.0035. Nowadays, the photon spectrum is 
confirmed to be a blackbody spectrum with a temperature given by I2T1 

T CMB = 2.725 ± 0.002 K (systematic, 95% c.l.) ± 7 /xK (la statistical) (98) 

In fact, this is the best blackbody spectrum ever measured, see Fig.|9j with spectral distortions below the 
level of 10 parts per million (ppm). 




Fig. 10: The Cosmic Microwave Background Spectrum seen by the DMR instrument on COBE. The top figure corresponds to 
the monopole, T = 2.725 ± 0.002 K. The middle figure shows the dipole, 5Ti = 3.372 ± 0.014 mK, and the lower figure 
shows the quadrupole and higher multipoles, 8T2 — 18 ± 2 fiK. The central region corresponds to foreground by the galaxy. 
From Ref. (23). 

Moreover, the differential microwave radiometer (DMR) instrument on COBE, with a resolution 
of about 7° in the sky, has also confirmed that it is an extraordinarily isotropic background. The devia- 
tions from isotropy, i.e. differences in the temperature of the blackbody spectrum measured in different 
directions in the sky, are of the order of 20 fiK on large scales, or one part in 10 5 , see Ref. 1251 . There 
is, in fact, a dipole anisotropy of one part in 10 3 , 5T\ = 3.372 ± 0.007 mK (95% c.l.), in the direction 
of the Virgo cluster, (I, b) = (264.14° ± 0.30, 48.26° ± 0.30) (95% c.l.). Under the assumption that a 
Doppler effect is responsible for the entire CMB dipole, the velocity of the Sun with respect to the CMB 



rest frame is v Q = 371 ± 0.5 km/s, see Ref. [21!|. 4 When subtracted, we are left with a whole spectrum 
of anisotropics in the higher multipoles (quadrupole, octupole, etc.), 6T2 = 18 ± 2 /jK (95% c.L), see 
Ref. G3 and Fig. GDI 

Soon after COBE, other groups quickly confirmed the detection of temperature anisotropics at 
around 30 /iK and above, at higher multipole numbers or smaller angular scales. As I shall discuss below, 
these anisotropics play a crucial role in the understanding of the origin of structure in the universe. 

2.69 Large-scale structure formation 

Although the isotropic microwave background indicates that the universe in the past was extraordinarily 
homogeneous, we know that the universe today is not exactly homogeneous: we observe galaxies, clus- 
ters and superclusters on large scales. These structures are expected to arise from very small primordial 
inhomogeneities that grow in time via gravitational instability, and that may have originated from tiny 
ripples in the metric, as matter fell into their troughs. Those ripples must have left some trace as temper- 
ature anisotropics in the microwave background, and indeed such anisotropics were finally discovered 
by the COBE satellite in 1992. The reason why they took so long to be discovered was that they appear 
as perturbations in temperature of only one part in 10 5 . 

While the predicted anisotropics have finally been seen in the CMB, not all kinds of matter and/or 
evolution of the universe can give rise to the structure we observe today. If we define the density contrast 
as |24) 



where p(a) = po a~ 3 is the average cosmic density, we need a theory that will grow a density contrast 
with amplitude 5 ~ 10~ 5 at the last scattering surface (z = 1100) up to density contrasts of the order of 
5 ~ 10 2 for galaxies at redshifts zCl, i.e. today. This is a necessary requirement for any consistent 
theory of structure formation [25 1. 

Furthermore, the anisotropics observed by the COBE satellite correspond to a small-amplitude 
scale-invariant primordial power spectrum of inhomogeneities 



where the brackets (•) represent integration over an ensemble of different universe realizations. These 
inhomogeneities are like waves in the space-time metric. When matter fell in the troughs of those waves, 
it created density perturbations that collapsed gravitationally to form galaxies and clusters of galaxies, 
with a spectrum that is also scale invariant. Such a type of spectrum was proposed in the early 1970s by 
Edward R. Harrison, and independently by the Russian cosmologist Yakov B. Zel'dovich, see Ref. [26 1, 
to explain the distribution of galaxies and clusters of galaxies on very large scales in our observable 
universe. 

Today various telescopes - like the Hubble Space Telescope, the twin Keck telescopes in Hawaii 
and the European Southern Observatory telescopes in Chile - are exploring the most distant regions of 
the universe and discovering the first galaxies at large distances. The furthest galaxies observed so far are 
at redshifts of z ~ 10 (at a distance of 13.7 billion light years from Earth), whose light was emitted when 
the universe had only about 3% of its present age. Only a few galaxies are known at those redshifts, but 
there are at present various catalogs like the Cf A and APM galaxy catalogs, and more recently the IRAS 
Point Source redshift Catalog, see Fig. ^2 and Las Campanas redshift surveys, that study the spatial 
distribution of hundreds of thousands of galaxies up to distances of a billion light years, or z < 0.1, 
or the 2 degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey (SDSS), 
which reach z < 0.5 and study millions of galaxies. These catalogs are telling us about the evolution 

4 COBE even determined the annual variation due to the Earth's motion around the Sun - the ultimate proof of Copernicus' 
hypothesis. 
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Fig. 11: The IRAS Point Source Catalog redshift survey contains some 15,000 galaxies, covering over 83% of the sky up to 
redshifts of z < 0.05. We show here the projection of the galaxy distribution in galactic coordinates. From Ref. 1271 . 

of clusters and superclusters of galaxies in the universe, and already put constraints on the theory of 
structure formation. From these observations one can infer that most galaxies formed at redshifts of the 
order of 2 — 6; clusters of galaxies formed at redshifts of order 1, and superclusters are forming now. 
That is, cosmic structure formed from the bottom up: from galaxies to clusters to superclusters, and not 
the other way around. This fundamental difference is an indication of the type of matter that gave rise to 
structure. 

We know from Big Bang nucleosynthesis that all the baryons in the universe cannot account for 
the observed amount of matter, so there must be some extra matter (dark since we don't see it) to account 
for its gravitational pull. Whether it is relativistic (hot) or non-relativistic (cold) could be inferred from 
observations: relativistic particles tend to diffuse from one concentration of matter to another, thus trans- 
ferring energy among them and preventing the growth of structure on small scales. This is excluded by 
observations, so we conclude that most of the matter responsible for structure formation must be cold. 
How much there is is a matter of debate at the moment. Some recent analyses suggest that there is not 
enough cold dark matter to reach the critical density required to make the universe flat. If we want to 
make sense of the present observations, we must conclude that some other form of energy permeates the 
universe. In order to resolve this issue, 2dFGRS and SDSS started taking data a few years ago. The first 
has already been completed, but the second one is still taking data up to redshifts z ~ 5 for quasars, over 
a large region of the sky. These important observations will help astronomers determine the nature of the 
dark matter and test the validity of the models of structure formation. 

Before COBE discovered the anisotropics of the microwave background there were serious doubts 
whether gravity alone could be responsible for the formation of the structure we observe in the universe 
today. It seemed that a new force was required to do the job. Fortunately, the anisotropics were found 
with the right amplitude for structure to be accounted for by gravitational collapse of primordial inho- 
mogeneities under the attraction of a large component of non-relativistic dark matter. Nowadays, the 
standard theory of structure formation is a cold dark matter model with a non vanishing cosmological 
constant in a spatially flat universe. Gravitational collapse amplifies the density contrast initially through 
linear growth and later on via non-linear collapse. In the process, overdense regions decouple from 
the Hubble expansion to become bound systems, which start attracting eachother to form larger bound 
structures. In fact, the largest structures, superclusters, have not yet gone non-linear. 

The primordial spectrum i 1001) is reprocessed by gravitational instability after the universe be- 
comes matter dominated and inhomogeneities can grow. Linear perturbation theory shows that the grow- 
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in the Einstein-de Sitter limit (u = p/p = 1/3 and 0, for radiation and matter, respectively). There are 
slight deviations for a 3> a cq , if 7^ 1 or S7a 7^ 0, but we will not be concerned with them here. 
The important observation is that, since the density contrast at last scattering is of order 5 ~ 1CT 5 , and 
the scale factor has grown since then only a factor z& ec ~ 10 3 , one would expect a density contrast 
today of order 5o ~ 1CT 2 . Instead, we observe structures like galaxies, where 5 ~ 10 2 . So how can 
this be possible? The microwave background shows anisotropics due to fluctuations in the baryonic 
matter component only (to which photons couple, electromagnetically). If there is an additional matter 
component that only couples through very weak interactions, fluctuations in that component could grow 
as soon as it decoupled from the plasma, well before photons decoupled from baryons. The reason why 
baryonic inhomogeneities cannot grow is because of photon pressure: as baryons collapse towards denser 
regions, radiation pressure eventually halts the contraction and sets up acoustic oscillations in the plasma 
that prevent the growth of perturbations, until photon decoupling. On the other hand, a weakly interacting 
cold dark matter component could start gravitational collapse much earlier, even before matter-radiation 
equality, and thus reach the density contrast amplitudes observed today. The resolution of this mismatch 
is one of the strongest arguments for the existence of a weakly interacting cold dark matter component 
of the universe. 
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Fig. 12: The power spectrum for cold dark matter (CDM), tilted cold dark matter (TCDM), hot dark matter (HDM), and mixed 
hot plus cold dark matter (MDM), normalized to COBE, for large-scale structure formation. From Ref. |28 1. 



How much dark matter there is in the universe can be deduced from the actual power spectrum (the 
Fourier transform of the two-point correlation function of density perturbations) of the observed large 
scale structure. One can decompose the density contrast in Fourier components, see Eq. d99t . This is 
very convenient since in linear perturbation theory individual Fourier components evolve independently. 
A comoving wavenumber k is said to "enter the horizon" when k = d~jl(a) = aH(a). If a certain 
perturbation, of wavelength A = k^ 1 < d//(a cq ), enters the horizon before matter-radiation equality, the 
fast radiation-driven expansion prevents dark-matter perturbations from collapsing. Since light can only 
cross regions that are smaller than the horizon, the suppression of growth due to radiation is restricted 
to scales smaller than the horizon, while large-scale perturbations remain unaffected. This is the reason 

5 The decaying modes go like 5(t) ~ t _1 , for all uj. 



why the horizon size at equality, Eq. d95t . sets an important scale for structure growth, 

k eq = ^(aeq) ^ 0.083 (n M h) h Mpc -1 . (102) 

The suppression factor can be easily computed from d 101l> as / sup = (a e nter/oeq) 2 = {k e q/k) 2 . In other 
words, the processed power spectrum P(k) will have the form: 




(103) 
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This is precisely the shape that large-scale galaxy catalogs are bound to test in the near future, see Fig. El 
Furthermore, since relativistic Hot Dark Matter (HDM) transfer energy between clumps of matter, they 
will wipe out small scale perturbations, and this should be seen as a distinctive signature in the matter 
power spectra of future galaxy catalogs. On the other hand, non-relativistic Cold Dark Matter (CDM) 
allow structure to form on all scales via gravitational collapse. The dark matter will then pull in the 
baryons, which will later shine and thus allow us to see the galaxies. 

Naturally, when baryons start to collapse onto dark matter potential wells, they will convert a large 
fraction of their potential energy into kinetic energy of protons and electrons, ionizing the medium. As a 
consequence, we expect to see a large fraction of those baryons constituting a hot ionized gas surrounding 
large clusters of galaxies. This is indeed what is observed, and confirms the general picture of structure 
formation. 



3. DETERMINATION OF COSMOLOGICAL PARAMETERS 

In this Section, I will restrict myself to those recent measurements of the cosmological parameters by 
means of standard cosmological techniques, together with a few instances of new results from recently 
applied techniques. We will see that a large host of observations are determining the cosmological 
parameters with some reliability of the order of 10%. However, the majority of these measurements are 
dominated by large systematic errors. Most of the recent work in observational cosmology has been 
the search for virtually systematic-free observables, like those obtained from the microwave background 
anisotropics, and discussed in Section 4.4. I will devote, however, this Section to the more 'classical' 
measurements of the following cosmological parameters: The rate of expansion Hq; the matter content 
Om; the cosmological constant S7a; the spatial curvature O^r, and the age of the universe to- 



3.1 The rate of expansion Hq 

Over most of last century the value of Ho has been a constant source of disagreement [29|. Around 
1929, Hubble measured the rate of expansion to be Ho = 500 kms _1 Mpc _1 , which implied an age of 
the universe of order to ~ 2 Gyr, in clear conflict with geology. Hubble's data was based on Cepheid 
standard candles that were incorrectly calibrated with those in the Large Magellanic Cloud. Later on, 
in 1954 Baade recalibrated the Cepheid distance and obtained a lower value, Hq = 250 kms _1 Mpc _1 , 
still in conflict with ratios of certain unstable isotopes. Finally, in 1958 Sandage realized that the bright- 
est stars in galaxies were ionized HII regions, and the Hubble rate dropped down to Ho = 60 kms -1 
Mpc -1 , still with large (factor of two) systematic errors. Fortunately, in the past 15 years there has 
been significant progress towards the determination of Ho, with systematic errors approaching the 10% 
level. These improvements come from two directions. First, technological, through the replacement of 
photographic plates (almost exclusively the source of data from the 1920s to 1980s) with charged couple 
devices (CCDs), i.e. solid state detectors with excellent flux sensitivity per pixel, which were previously 
used successfully in particle physics detectors. Second, by the refinement of existing methods for mea- 
suring extragalactic distances (e.g. parallax, Cepheids, supernovae, etc.). Finally, with the development 
of completely new methods to determine Hq, which fall into totally independent and very broad cate- 
gories: a) Gravitational lensing; b) Sunyaev-Zel'dovich effect; c) Extragalactic distance scale, mainly 



Cepheid variability and type la Supernovae; d) Microwave background anisotropics. I will review here 
the first three, and leave the last method for Section 4.4, since it involves knowledge about the primordial 
spectrum of inhomogeneities. 

3.11 Gravitational lensing 

Imagine a quasi-stellar object (QSO) at large redshift (z 3> 1) whose light is lensed by an intervening 
galaxy at redshift z ~ 1 and arrives to an observer at z = 0. There will be at least two different 
images of the same background variable point source. The arrival times of photons from two different 
gravitationally lensed images of the quasar depend on the different path lengths and the gravitational 
potential traversed. Therefore, a measurement of the time delay and the angular separation of the different 
images of a variable quasar can be used to determine Hq with great accuracy. This method, proposed in 
1964 by Refsdael [30], offers tremendous potential because it can be applied at great distances and it is 
based on very solid physical principles l3Tll . 

Unfortunately, there are very few systems with both a favourable geometry (i.e. a known mass 
distribution of the intervening galaxy) and a variable background source with a measurable time delay. 
That is the reason why it has taken so much time since the original proposal for the first results to come 
out. Fortunately, there are now very powerful telescopes that can be used for these purposes. The best 
candidate to-date is the QSO 0957 + 561, observed with the 10m Keck telescope, for which there is a 
model of the lensing mass distribution that is consistent with the measured velocity dispersion. Assuming 
a flat space with Qm = 0.25, one can determine 1321 

H = 72 ± 7 (la statistical) ± 15% (systematic) kms^Mpc -1 . (104) 

The main source of systematic error is the degeneracy between the mass distribution of the lens and 
the value of Hq. Knowledge of the velocity dispersion within the lens as a function of position helps 
constrain the mass distribution, but those measurements are very difficult and, in the case of lensing by 
a cluster of galaxies, the dark matter distribution in those systems is usually unknown, associated with a 
complicated cluster potential. Nevertheless, the method is just starting to give promising results and, in 
the near future, with the recent discovery of several systems with optimum properties, the prospects for 
measuring Hq and lowering its uncertainty with this technique are excellent. 

3.12 Sunyaev-Zel'dovich effect 

As discussed in the previous Section, the gravitational collapse of baryons onto the potential wells gen- 
erated by dark matter gave rise to the reionization of the plasma, generating an X-ray halo around rich 
clusters of galaxies, see Fig. ^] The inverse-Compton scattering of microwave background photons 
off the hot electrons in the X-ray gas results in a measurable distortion of the blackbody spectrum 
of the microwave background, known as the Sunyaev-Zel'dovich (SZ) effect. Since photons acquire 
extra energy from the X-ray electrons, we expect a shift towards higher frequencies of the spectrum, 
(Au/v) ~ (k-QT Z3S /m e c 2 ) ~ 10~ 2 . This corresponds to a decrement of the microwave background tem- 
perature at low frequencies (Ray leigh- Jeans region) and an increment at high frequencies, see Ref. 1331 . 

Measuring the spatial distribution of the SZ effect (3 K spectrum), together with a high resolution 
X-ray map (10 8 K spectrum) of the cluster, one can determine the density and temperature distribution 
of the hot gas. Since the X-ray flux is distance-dependent (J 7 = C/Aird 2 L ), while the SZ decrement is 
not (because the energy of the CMB photons increases as we go back in redshift, u = uq(1 + z), and 
exactly compensates the redshift in energy of the photons that reach us), one can determine from there 
the distance to the cluster, and thus the Hubble rate Hq. 

The advantages of this method are that it can be applied to large distances and it is based on 
clear physical principles. The main systematics come from possible dumpiness of the gas (which would 



Fig. 13: The Coma cluster of galaxies, seen here in an optical image (left) and an X-ray image (right), taken by the recently 
launched Chandra X-ray Observatory. From Ref. [34|. 

reduce Ho), projection effects (if the clusters are prolate, Ho could be larger), the assumption of hy- 
drostatic equilibrium of the X-ray gas, details of models for the gas and electron densities, and possible 
contaminations from point sources. Present measurements give the value [33 1 

H = 60 ± 10 (la statistical) ± 20% (systematic) kms^Mpc -1 , (105) 

compatible with other determinations. A great advantage of this completely new and independent method 
is that nowadays more and more clusters are observed in the X-ray, and soon we will have high-resolution 
2D maps of the SZ decrement from several balloon flights, as well as from future microwave background 
satellites, together with precise X-ray maps and spectra from the Chandra X-ray observatory recently 
launched by NASA, as well as from the European X-ray satellite XMM launched a few months ago by 
ESA, which will deliver orders of magnitude better resolution than the existing Einstein X-ray satellite. 

3.13 Cepheid variability 

Cepheids are low-mass variable stars with a period-luminosity relation based on the helium ionization 
cycles inside the star, as it contracts and expands. This time variability can be measured, and the star's 
absolute luminosity determined from the calibrated relationship. From the observed flux one can then 
deduce the luminosity distance, see Eq. d28l . and thus the Hubble rate Hq. The Hubble Space Telescope 
(HST) was launched by NASA in 1990 (and repaired in 1993) with the specific project of calibrating the 
extragalactic distance scale and thus determining the Hubble rate with 10% accuracy. The most recent 
results from HST are the following [35 1 

Ho = 71 ± 4 (random) ± 7 (systematic) kms _1 Mpc~ 1 . (106) 

The main source of systematic error is the distance to the Large Magellanic Cloud, which provides the 
fiducial comparison for Cepheids in more distant galaxies. Other systematic uncertainties that affect the 
value of Ho are the internal extinction correction method used, a possible metallicity dependence of the 
Cepheid period-luminosity relation and cluster population incompleteness bias, for a set of 2 1 galaxies 
within 25 Mpc, and 23 clusters within z % 0.03. 

With better telescopes already taking data, like the Very Large Telescope (VLT) interferometer 
of the European Southern Observatory (ESO) in the Chilean Atacama desert, with 8 synchronized tele- 
scopes, and others coming up soon, like the Next Generation Space Telescope (NGST) proposed by 
NASA for 2008, and the Gran TeCan of the European Northern Observatory in the Canary Islands, for 
2010, it is expected that much better resolution and therefore accuracy can be obtained for the determi- 
nation of Hq- 



3.2 Dark Matter 

In the 1920s Hubble realized that the so called nebulae were actually distant galaxies very similar to our 
own. Soon afterwards, in 1933, Zwicky found dynamical evidence that there is possibly ten to a hundred 
times more mass in the Coma cluster than contributed by the luminous matter in galaxies [ 36 1. However, 
it was not until the 1970s that the existence of dark matter began to be taken more seriously. At that time 
there was evidence that rotation curves of galaxies did not fall off with radius and that the dynamical 
mass was increasing with scale from that of individual galaxies up to clusters of galaxies. Since then, 
new possible extra sources to the matter content of the universe have been accumulating: 
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The empirical route to the determination of Qm is nowadays one of the most diversified of all 
cosmological parameters. The matter content of the universe can be deduced from the mass-to-light ratio 
of various objects in the universe; from the rotation curves of galaxies; from microlensing and the direct 
search of Massive Compact Halo Objects (MACHOs); from the cluster velocity dispersion with the use 
of the Virial theorem; from the baryon fraction in the X-ray gas of clusters; from weak gravitational 
lensing; from the observed matter distribution of the universe via its power spectrum; from the cluster 
abundance and its evolution; from direct detection of massive neutrinos at SuperKamiokande; from direct 
detection of Weakly Interacting Massive Particles (WIMPs) at CDMS, DAMA or UKDMC, and finally 
from microwave background anisotropics. I will review here just a few of them. 

3.21 Rotation curves of spiral galaxies 

The flat rotation curves of spiral galaxies provide the most direct evidence for the existence of large 
amounts of dark matter. Spiral galaxies consist of a central bulge and a very thin disk, stabilized against 
gravitational collapse by angular momentum conservation, and surrounded by an approximately spher- 
ical halo of dark matter. One can measure the orbital velocities of objects orbiting around the disk as a 
function of radius from the Doppler shifts of their spectral lines. 

The rotation curve of the Andromeda galaxy was first measured by Babcock in 1938, from the 
stars in the disk. Later it became possible to measure galactic rotation curves far out into the disk, and 
a trend was found [37|. The orbital velocity rose linearly from the center outward until it reached a 
typical value of 200 km/s, and then remained flat out to the largest measured radii. This was completely 
unexpected since the observed surface luminosity of the disk falls off exponentially with radius 1371 . 
I(r) = Jq exp(— r/ro). Therefore, one would expect that most of the galactic mass is concentrated 
within a few disk lengths rr>, such that the rotation velocity is determined as in a Keplerian orbit, v mt = 
(GM/r) 1 / 2 oc r -1 / 2 . No such behaviour is observed. In fact, the most convincing observations come 
from radio emission (from the 21 cm line) of neutral hydrogen in the disk, which has been measured to 
much larger galactic radii than optical tracers. A typical case is that of the spiral galaxy NGC 6503, where 
rr> = 1.73 kpc, while the furthest measured hydrogen line is at r = 22.22 kpc, about 13 disk lengths 
away. Nowadays, thousands of galactic rotation curves are known, see Fig. 14, and all suggest the 
existence of about ten times more mass in the halos of spiral galaxies than in the stars of the disk. Recent 
numerical simulations of galaxy formation in a CDM cosmology [38] suggest that galaxies probably 
formed by the infall of material in an overdense region of the universe that had decoupled from the 
overall expansion. 

The dark matter is supposed to undergo violent relaxation and create a virialized system, i.e. 
in hydrostatic equilibrium. This picture has led to a simple model of dark-matter halos as isothermal 
spheres, with density profile p(r) = p c /{r1 + r 2 ), where r c is a core radius and p c = v^/AttG, with 




Fig. 14: The rotation curves of several hundred galaxies. Upper panel: As a function of their radii in kpc. Middle panel: The 
central 5 kpc. Lower panel: As a function of scale radius. 



Voo equal to the plateau value of the flat rotation curve. This model is consistent with the universal 
rotation curves seen in Fig. 6. At large radii the dark matter distribution leads to a flat rotation curve. 
The question is for how long. In dense galaxy clusters one expects the galactic halos to overlap and 
form a continuum, and therefore the rotation curves should remain flat from one galaxy to another. 
However, in field galaxies, far from clusters, one can study the rotation velocities of substructures (like 
satellite dwarf galaxies) around a given galaxy, and determine whether they fall off at sufficiently large 
distances according to Kepler's law, as one would expect, once the edges of the dark matter halo have 
been reached. These observations are rather difficult because of uncertainties in distinguishing between 
true satellites and interlopers. Recently, a group from the Sloan Digital Sky Survey Collaboration claim 
that they have seen the edges of the dark matter halos around field galaxies by confirming the fall-off 
at large distances of their rotation curves [39]. These results, if corroborated by further analysis, would 
constitute a tremendous support to the idea of dark matter as a fluid surrounding galaxies and clusters, 
while at the same time eliminates the need for modifications of Newtonian of even Einstenian gravity at 
the scales of galaxies, to account for the flat rotation curves. 

That's fine, but how much dark matter is there at the galactic scale? Adding up all the matter in 
galactic halos up to a maximum radii, one finds 

ft halo ~ 10 n lum > 0.03 - 0.05 . (Ill) 

Of course, it would be extraordinary if we could confirm, through direct detection, the existence of dark 
matter in our own galaxy. For that purpose, one should measure its rotation curve, which is much more 



difficult because of obscuration by dust in the disk, as well as problems with the determination of reliable 
galactocentric distances for the tracers. Nevertheless, the rotation curve of the Milky Way has been 
measured and conforms to the usual picture, with a plateau value of the rotation velocity of 220 km/s. 
For dark matter searches, the crucial quantity is the dark matter density in the solar neighbourhood, which 
turns out to be (within a factor of two uncertainty depending on the halo model) poM = 0-3 GeV/cm 3 . 
We will come back to direct searched of dark matter in a later subsection. 

3.22 Baryon fraction in clusters 

Since large clusters of galaxies form through gravitational collapse, they scoop up mass over a large 
volume of space, and therefore the ratio of baryons over the total matter in the cluster should be rep- 
resentative of the entire universe, at least within a 20% systematic error. Since the 1960s, when X-ray 
telescopes became available, it is known that galaxy clusters are the most powerful X-ray sources in the 
sky BOl . The emission extends over the whole cluster and reveals the existence of a hot plasma with 
temperature T ~ 10 7 — 10 8 K, where X-rays are produced by electron bremsstrahlung. Assuming the 
gas to be in hydrostatic equilibrium and applying the virial theorem one can estimate the total mass in 
the cluster, giving general agreement (within a factor of 2) with the virial mass estimates. From these 
estimates one can calculate the baryon fraction of clusters 

f B h 3/2 = 0.08 => ^£ « 0.14, for ft = 0.70. (112) 

Since Q\ um ~ 0.002 — 0.006, the previous expression suggests that clusters contain far more baryonic 
matter in the form of hot gas than in the form of stars in galaxies. Assuming this fraction to be repre- 
sentative of the entire universe, and using the Big Bang nucleosynthesis value of Qb = 0.04 ± 0.01, for 
h = 0.7, we find 

n M = 0.3 ± 0.1 (statistical) ± 20% (systematic) . (113) 

This value is consistent with previous determinations of Qm- If some baryons are ejected from the cluster 
during gravitational collapse, or some are actually bound in nonluminous objects like planets, then the 
actual value of Qm is smaller than this estimate. 

3.23 Weak gravitational lensing 

Since the mid 1980s, deep surveys with powerful telescopes have observed huge arc-like features in 
galaxy clusters. The spectroscopic analysis showed that the cluster and the giant arcs were at very differ- 
ent redshifts. The usual interpretation is that the arc is the image of a distant background galaxy which 
is in the same line of sight as the cluster so that it appears distorted and magnified by the gravitational 
lens effect: the giant arcs are essentially partial Einstein rings. From a systematic study of the clus- 
ter mass distribution one can reconstruct the shear field responsible for the gravitational distortion 1411 . 
This analysis shows that there are large amounts of dark matter in the clusters, in rough agreement with 
the virial mass estimates, although the lensing masses tend to be systematically larger. At present, the 
estimates indicate Qm = 0.2 — 0.3 on scales % 6 h^ 1 Mpc. 

3.24 Large scale structure formation and the matter power spectrum 

Although the isotropic microwave background indicates that the universe in the past was extraordinarily 
homogeneous, we know that the universe today is far from homogeneous: we observe galaxies, clusters 
and superclusters on large scales. These structures are expected to arise from very small primordial inho- 
mogeneities that grow in time via gravitational instability, and that may have originated from tiny ripples 
in the metric, as matter fell into their troughs. Those ripples must have left some trace as temperature 
anisotropics in the microwave background, and indeed such anisotropics were finally discovered by the 





2dF Galaxy Redshift Survey 




Fig. 15: The 2 degree Field Galaxy Redshift Survey contains some 250,000 galaxies, covering a large fraction of the sky up to 
redshifts of z < 0.25. From Ref. |42|. 



COBE satellite in 1992. However, not all kinds of matter and/or evolution of the universe can give rise 
to the structure we observe today. If we define the density contrast as 



where p(a) = po a~ 3 is the average cosmic density, we need a theory that will grow a density contrast 
with amplitude 5 ~ 1CP 5 at the last scattering surface (z = 1100) up to density contrasts of the order of 
5 ~ 10 2 for galaxies at redshifts i.e. today. This is a necessary requirement for any consistent 

theory of structure formation. 

Furthermore, the anisotropics observed by the COBE satellite correspond to a small-amplitude 
scale-invariant primordial power spectrum of inhomogeneities 



These inhomogeneities are like waves in the space-time metric. When matter fell in the troughs of those 
waves, it created density perturbations that collapsed gravitationally to form galaxies and clusters of 
galaxies, with a spectrum that is also scale invariant. Such a type of spectrum was proposed in the early 
1970s by Edward R. Harrison, and independently by the Russian cosmologist Yakov B. Zel'dovich [26 1, 
to explain the distribution of galaxies and clusters of galaxies on very large scales in our observable 
universe, see Fig. 15. 

Since the primordial spectrum is very approximately represented by a scale-invariant Gaussian 
random field, the best way to present the results of structure formation is by working with the 2-point 
correlation function in Fourier space, the so-called power spectrum. If the reprocessed spectrum of in- 
homogeneities remains Gaussian, the power spectrum is all we need to describe the galaxy distribution. 
Non-Gaussian effects are expected to arise from the non-linear gravitational collapse of structure, and 
may be important at small scales. The power spectrum measures the degree of inhomogeneity in the 
mass distribution on different scales, see Fig. 16. It depends upon a few basic ingredientes: a) the pri- 
mordial spectrum of inhomogeneities, whether they are Gaussian or non-Gaussian, whether adiabatic 
(perturbations in the energy density) or isocurvature (perturbations in the entropy density), whether the 
primordial spectrum has tilt (deviations from scale-invariance), etc.; b) the recent creation of inhomo- 
geneities, whether cosmic strings or some other topological defect from an early phase transition are 
responsible for the formation of structure today; and c) the cosmic evolution of the inhomogeneity, 
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Fig. 16: The measured power spectrum P(k) as a function of wavenumber k. From observations of the Sloan Digital Sky 
Survey, CMB anisotropics, cluster abundance, gravitational lensing and Lyman-a forest. From Ref. |43|. 



whether the universe has been dominated by cold or hot dark matter or by a cosmological constant since 
the beginning of structure formation, and also depending on the rate of expansion of the universe. 

The working tools used for the comparison between the observed power spectrum and the pre- 
dicted one are very precise N-body numerical simulations and theoretical models that predict the shape 
but not the amplitude of the present power spectrum. Even though a large amount of work has gone 
into those analyses, we still have large uncertainties about the nature and amount of matter necessary for 
structure formation. A model that has become a working paradigm is a flat cold dark matter model with 
a cosmological constant and &m ~ 0-3- This model is now been confronted with the recent very precise 
measurements from 2dFGRS E3 and SDSS 1431 . 

3.25 The new redshift catalogs, 2dF and Sloan Digital Sky Survey 

Our view of the large-scale distribution of luminous objects in the universe has changed dramatically 
during the last 25 years: from the simple pre- 1975 picture of a distribution of field and cluster galax- 
ies, to the discovery of the first single superstructures and voids, to the most recent results showing an 
almost regular web-like network of interconnected clusters, filaments and walls, separating huge nearly 
empty volumes. The increased efficiency of redshift surveys, made possible by the development of spec- 
trographs and - specially in the last decade - by an enormous increase in multiplexing gain (i.e. the 
ability to collect spectra of several galaxies at once, thanks to fibre-optic spectrographs), has allowed 
us not only to do cartography of the nearby universe, but also to statistically characterize some of its 
properties. At the same time, advances in theoretical modeling of the development of structure, with 
large high-resolution gravitational simulations coupled to a deeper yet limited understanding of how to 
form galaxies within the dark matter halos, have provided a more realistic connection of the models to 
the observable quantities. Despite the large uncertainties that still exist, this has transformed the study of 
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Fig. 17: The observed cosmic matter components as functions of the Hubble expansion parameter. The luminous matter 
component is given by 0.002 < fii um < 0.006; the galactic halo component is the horizontal band, 0.03 < fihaio < 0.05, 
crossing the baryonic component from BBN, Qb h 2 = 0.0244 ± 0.0024; and the dynamical mass component from large scale 
structure analysis is given by Qm = 0.3 ± 0.1. Note that in the range Ho = 70 ± 7 km/s/Mpc, there are three dark matter 
problems, see the text. From Ref. |44|. 



cosmology and large-scale structure into a truly quantitative science, where theory and observations can 
progress together. 

3.26 Summary of the matter content 

We can summarize the present situation with Fig. 17, for Q,m as a function of Ho. There are four bands, 
the luminous matter £l\ um ; the baryon content from BBN; the galactic halo component Slhaio 5 and 
the dynamical mass from clusters, £Im- From this figure it is clear that there are in fact three dark matter 
problems: The first one is where are 90% of the baryons? Between the fraction predicted by BBN and 
that seen in stars and diffuse gas there is a huge fraction which is in the form of dark baryons. They could 
be in small clumps of hydrogen that have not started thermonuclear reactions and perhaps constitute the 
dark matter of spiral galaxies' halos. Note that although Qb and £\aio coincide at Hq ~ 70 km/s/Mpc, 
this could be just a coincidence. The second problem is what constitutes 90% of matter, from BBN 
baryons to the mass inferred from cluster dynamics? This is the standard dark matter problem and could 
be solved in the future by direct detection of a weakly interacting massive particle in the laboratory. And 
finally, since we know from observations of the CMB that the universe is flat, the rest, up to Q,q = 1, 
must be a diffuse vacuum energy, which affects the very large scales and late times, and seems to be 
responsible for the present acceleration of the universe, see Section 3. Nowadays, multiple observations 
seem to converge towards a common determination of Qm = 0.25 ± 0.08 (95% c.l.), see Fig. 18. 

3.27 Massive neutrinos 

One of the 'usual suspects' when addressing the problem of dark matter are neutrinos. They are the only 
candidates known to exist. If neutrinos have a mass, could they constitute the missing matter? We know 
from the Big Bang theory, see Section 2.6.5, that there is a cosmic neutrino background at a temperature 
of approximately 2K. This allows one to compute the present number density in the form of neutrinos, 
which turns out to be, for massless neutrinos, n v (T v ) = n 7 (T 7 ) = 112 cm -3 , per species of neutrino. 
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Fig. 18: Different determinations of S1m as a function of distance, from various sources: 1. peculiar velocities; 2. weak 
gravitacional lensing; 3. shear autocorrelation function; 4. local group of galaxies; 5. baryon mass fraction; 6. cluster mass 
function; 7. virgocentric flow; 8. mean relative velocities; 9. redshift space distortions; 10. mass power spectrum; 11. integrated 
Sachs- Wolfe effect; 12. angular diameter distance: SNe; 13. cluster baryon fraction. While a few years ago the dispersion 
among observed values was huge and strongly dependent on scale, at present the observed value of the matter density parameter 
falls well within a narrow range, Qm = 0.25 ± 0.07 (95% c.l.) and is essentially independent on scale, from 100 kpc to 5000 
Mpc. Adapted from Ref. l45l . 
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If neutrinos have mass, as recent experiments seem to suggest, see Fig. 19, the cosmic energy density 

IT n 7 



in massive neutrinos would be p v = ^ n^m^ = £■ n 7 J2 m w an d therefore its contribution today, 



n v h 2 = ^%. (116) 
93.2 eV 

The discussion in the previous Sections suggest that < 0.4, and thus, for any of the three families 
of neutrinos, m v < 40 eV. Note that this limit improves by six orders of magnitude the present bound 
on the tau-neutrino mass fl9"l . Supposing that the missing mass in non-baryonic cold dark matter arises 
from a single particle dark matter (PDM) component, its contribution to the critical density is bounded 
by 0.05 < Qp DM h 2 < 0.4, see Fig. 17. 

I will now go through the various logical arguments that exclude neutrinos as the dominant compo- 
nent of the missing dark matter in the universe. Is it possible that neutrinos with a mass 4 eV < m u < 40 
eV be the non-baryonic PDM component? For instance, could massive neutrinos constitute the dark 
matter halos of galaxies? For neutrinos to be gravitationally bound to galaxies it is necessary that their 
velocity be less that the escape velocity v csc , and thus their maximum momentum is p max = m u v esc . 
How many neutrinos can be packed in the halo of a galaxy? Due to the Pauli exclusion principle, 
the maximum number density is given by that of a completely degenerate Fermi gas with momen- 
tum pf = p max , i.e. n max = £>m ax /37r 2 . Therefore, the maximum local density in dark matter 
neutrinos is p ma , x = n ma _ x m u = mf,v^ sc /37r 2 , which must be greater than the typical halo density 
Phalo = 0.3 GeVcm -3 . For a typical spiral galaxy, this constraint, known as the Tremaine-Gunn limit, 
gives m v > 40 eV, see Ref. 1471 . However, this mass, even for a single species, say the tau-neutrino, 
gives a value for VL v h 2 = 0.5, which is far too high for structure formation. Neutrinos of such a low 
mass would constitute a relativistic hot dark matter component, which would wash-out structure below 
the supercluster scale, against evidence from present observations, see Fig. 19. Furthermore, apply- 



6 For a review on Neutrino properties, see Gonzalez-Garcfa's lectures on these Proceedings. 
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Fig. 19: The neutrino parameter space, mixing angle against Am 2 , including the results from the different solar and atmospheric 
neutrino oscillation experiments. Note the threshold of cosmologically important masses, cosmologically detectable neutrinos 
(by CMB and LSS observations), and cosmologically excluded range of masses. Adapted from Refs. |46 1 and |91 1. 



ing the same phase-space argument to the neutrinos as dark matter in the halo of dwarf galaxies gives 
m v > 100 eV, beyond closure density (U16t . We must conclude that the simple idea that light neutrinos 
could constitute the particle dark matter on all scales is ruled out. They could, however, still play a role 
as a sub-dominant hot dark matter component in a flat CDM model. In that case, a neutrino mass of order 
1 eV is not cosmological excluded, see Fig. 19. 

Another possibility is that neutrinos have a large mass, of order a few GeV. In that case, their num- 
ber density at decoupling, see Section 2.5.1, is suppressed by a Boltzmann factor, ~ exp(— m v /T^ ec ). 
For masses m v > T^ cc ~ 0.8 MeV, the present energy density has to be computed as a solution 
of the corresponding Boltzmann equation. Apart from a logarithmic correction, one finds £l v h 2 ~ 
0.1(10 GeV/m u ) 2 for Majorana neutrinos and slightly smaller for Dirac neutrinos. In either case, neu- 
trinos could be the dark matter only if their mass was a few GeV. Laboratory limits for v T of around 1 8 
MeV [19], and much more stringent ones for and u e , exclude the known light neutrinos. However, 
there is always the possibility of a fourth unknown heavy and stable (perhaps sterile) neutrino. If it 
couples to the Z boson and has a mass below 45 GeV for Dirac neutrinos (39.5 GeV for Majorana neu- 
trinos), then it is ruled out by measurements at LEP of the invisible width of the Z. There are two logical 
alternatives, either it is a sterile neutrino (it does not couple to the Z), or it does couple but has a larger 
mass. In the case of a Majorana neutrino (its own antiparticle), their abundance, for this mass range, 
is too small for being cosmologically relevant, Q v h 2 < 0.005. If it were a Dirac neutrino there could 
be a lepton asymmetry, which may provide a higher abundance (similar to the case of baryogenesis). 
However, neutrinos scatter on nucleons via the weak axial-vector current (spin-dependent) interaction. 
For the small momentum transfers imparted by galactic WIMPs, such collisions are essentially coherent 
over an entire nucleus, leading to an enhancement of the effective cross section. The relatively large 
detection rate in this case allowes one to exclude fourth-generation Dirac neutrinos for the galactic dark 
matter [48 ]. Anyway, it would be very implausible to have such a massive neutrino today, since it would 
have to be stable, with a life-time greater than the age of the universe, and there is no theoretical reason 



to expect a massive sterile neutrino that does not oscillate into the other neutrinos. 

Of course, the definitive test to the possible contribution of neutrinos to the overall density of 
the universe would be to measure directly their mass in laboratory experiments. There are at present 
two types of experiments: neutrino oscillation experiments, which measure only differences in squared 
masses, and direct mass-searches experiments, like the tritium /3-spectrum and the neutrinoless double-/? 
decay experiments, which measure directly the mass of the electron neutrino. The former experiments 
give a bound m Ue ;$ 2.3 eV (95% c.l.) [49), while the latter claim [ 50 1 they have a positive evidence 
for a Majorana neutrino of mass m„ = 0.05 — 0.89 eV (95% c.l.), although this result still awaits 
confirmation by other experiments. Neutrinos with such a mass could very well constitute the HDM 
component of the universe, Ohdm ~ 0.15. The oscillation experiments give a range of possibilities 
for Am 2 = 0.3 — 3 eV 2 from LSND (not yet confirmed by Miniboone), to the atmospheric neutrino 
oscillations from Super Kamiokande (Am 2 ~ 2.2 ± 0.5 x 10~ 3 eV 2 , tan 2 6 = 1.0 ± 0.3) and the 
solar neutrino oscillations from KamLAND and the Sudbury Neutrino Observatory (Am 2 ~ 8.2 ± 
0.3 x 10~ 5 eV 2 , tan 2 6> = 0.39 ± 0.05), see Ref. 146). Only the first two possibilities would be 
cosmologically relevant, see Fig. 19. Thanks to recent observations by WMAP, 2dFGRS and SDSS, we 
can put stringent limits on the absolute scale of neutrino masses, see below (Section 3.4). 



3.28 Weakly Interacting Massive Particles 

Unless we drastically change the theory of gravity on large scales, baryons cannot make up the bulk 
of the dark matter. Massive neutrinos are the only alternative among the known particles, but they are 
essentially ruled out as a universal dark matter candidate, even if they may play a subdominant role as 
a hot dark matter component. There remains the mystery of what is the physical nature of the dominant 
cold dark matter component. Something like a heavy stable neutrino, a generic Weakly Interacting 
Massive Particle (WIMP), could be a reasonable candidate because its present abundance could fall 
within the expected range, 
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Here v IC \ is the relative velocity of the two incoming dark matter particles and the brackets (•) denote a 
thermal average at the freeze-out temperature, Tf ~ mprjivi/20, when the dark matter particles go out 
of equilibrium with radiation. The value of (cr ann ?; re i) needed for Opdm ~ 1 is remarkably close to 
what one would expect for a WIMP with a mass mpoM = 100 GeV, ((Xann^rci) ~ a 2 /8?TmpDM ~ 
3 x 10 -27 cm 3 s _1 . We still do not know whether this is just a coincidence or an important hint on the 
nature of dark matter. 

There are a few theoretical candidates for WIMPs, like the neutralino, coming from supersymme- 
tric extensions of the standard model of particle physics, 7 but at present there is no empirical evidence that 
such extensions are indeed realized in nature. In fact, the non-observation of supersymmetric particles 
at current accelerators places stringent limits on the neutralino mass and interaction cross section [52|. 
If WIMPs constitute the dominant component of the halo of our galaxy, it is expected that some may 
cross the Earth at a reasonable rate to be detected. The direct experimental search for them rely on 
elastic WIMP collisions with the nuclei of a suitable target. Dark matter WIMPs move at a typical 
galactic "virial" velocity of around 200 — 300 km/s, depending on the model. If their mass is in the 
range 10 — 100 GeV, the recoil energy of the nuclei in the elastic collision would be of order 10 keV. 
Therefore, one should be able to identify such energy depositions in a macroscopic sample of the target. 
There are at present three different methods: First, one could search for scintillation light in Nal crystals 
or in liquid xenon; second, search for an ionization signal in a semiconductor, typically a very pure 
germanium crystal; and third, use a cryogenic detector at 10 mK and search for a measurable temperature 



7 For a review of Supersymmetry (SUSY), see Kazakov's contribution to these Proceedings. 
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Fig. 20: The annual-modulation signal accumulated over 7 years is consistent with a neutralino of mass of m x 
and a proton cross section of £a p — 7.0 ~JZ? t x 10~ 6 pb, according to DAMA. From Ref. 1511 . 
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increase of the sample. The main problem with such a type of experiment is the low expected signal rate, 
with a typical number below 1 event/kg/day. To reduce natural radioactive contamination one must 
use extremely pure substances, and to reduce the background caused by cosmic rays requires that these 
experiments be located deeply underground. 

The best limits on WIMP scattering cross sections come from some germanium experiments, like 
the Criogenic Dark Matter Search (CDMS) collaboration at Stanford and the Soudan mine 1531 . as well 
as from the Nal scintillation detectors of the UK dark matter collaboration (UKDMC) in the Boulby salt 
mine in England [54], and the DAMA experiment in the Gran Sasso laboratory in Italy IBTll . Current 
experiments already touch the parameter space expected from supersymmetric particles, see Fig. 21, 
and therefore there is a chance that they actually discover the nature of the missing dark matter. The 
problem, of course, is to attribute a tentative signal unambiguously to galactic WIMPs rather than to 
some unidentified radioactive background. 

One specific signature is the annual modulation which arises as the Earth moves around the Sun. 8 
Therefore, the net speed of the Earth relative to the galactic dark matter halo varies, causing a modulation 
of the expected counting rate. The DAMA/Nal experiment has actually reported such a modulation 
signal, from the combined analysis of their 7-year data, see Fig. 20 and Ref. [51 1, which provides a 
confidence level of 99.6% for a neutralino mass of m x = 52 lg° GeV and a proton cross section of 
&p = 7.2 t° .t x 10" 6 pb, where £ = p x /0.3 GeV cm 3 is the local neutralino energy density in 
units of the galactic halo density. There has been no confirmation yet of this result from other dark 
matter search groups. In fact, the CDMS collaboration claims an exclusion of the DAMA region at 
the 3 sigma level, see Fig. 21. Hopefully in the near future we will have much better sensitivity at 
low masses from the Cryogenic Rare Event Search with Superconducting Thermometers (CRESST) 
experiment at Gran Sasso. The CRESST experiment [55] uses sapphire crystals as targets and a new 
method to simultaneously measure the phonons and the scintillating light from particle interactions inside 



The time scale of the Sun's orbit around the center of the galaxy is too large to be relevant in the analysis. 
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Fig. 21: Exclusion range for the spin-independent WIMP scattering cross section per nucleon from the Nal experiments and 
the Ge detectors. The blue lines come from the CDMS experiment, which exclude the DAMA region at more than 3 sigma. 
Also shown in yellow and red is the range of expected counting rates for neutralinos in the MSSM. From Ref. 1 53 . 



the crystal, which allows excellent background discrimination. Very recently there has been also the 
proposal of a completely new method based on a Superheated Droplet Detector (SDD), which claims to 
have already a similar sensitivity as the more standard methods described above, see Ref. J56j. 

There exist other indirect methods to search for galactic WIMPs [57]. Such particles could self- 
annihilate at a certain rate in the galactic halo, producing a potentially detectable background of high 
energy photons or antiprotons. The absence of such a background in both gamma ray satellites and 
the Alpha Matter Spectrometer [58 1 imposes bounds on their density in the halo. Alternatively, WIMPs 
traversing the solar system may interact with the matter that makes up the Earth or the Sun so that a small 
fraction of them will lose energy and be trapped in their cores, building up over the age of the universe. 
Their annihilation in the core would thus produce high energy neutrinos from the center of the Earth or 
from the Sun which are detectable by neutrino telescopes. In fact, Super Kamiokande already covers a 
large part of SUSY parameter space. In other words, neutrino telescopes are already competitive with 
direct search experiments. In particular, the AMANDA experiment at the South Pole [59], which has 
approximately 10 3 Cherenkov detectors several km deep in very clear ice, over a volume ~ 1 km 3 , is 
competitive with the best direct searches proposed. The advantages of AMANDA are also directional, 
since the arrays of Cherenkov detectors will allow one to reconstruct the neutrino trajectory and thus its 
source, whether it comes from the Earth or the Sun. AMANDA recently reported the detection of TeV 
neutrinos [59|. 

3.3 The age of the universe to 

The universe must be older than the oldest objects it contains. Those are believed to be the stars in the 
oldest clusters in the Milky Way, globular clusters. The most reliable ages come from the application 
of theoretical models of stellar evolution to observations of old stars in globular clusters. For about 30 
years, the ages of globular clusters have remained reasonable stable, at about 15 Gyr [ 60 1 . However, 
recently these ages have been revised downward [61 1. 

During the 1980s and 1990s, the globular cluster age estimates have improved as both new obser- 



vations have been made with CCDs, and since refinements to stellar evolution models, including opaci- 
ties, consideration of mixing, and different chemical abundances have been incorporated [62]. From the 
theory side, uncertainties in globular cluster ages come from uncertainties in convection models, opac- 
ities, and nuclear reaction rates. From the observational side, uncertainties arise due to corrections for 
dust and chemical composition. However, the dominant source of systematic errors in the globular clus- 
ter age is the uncertainty in the cluster distances. Fortunately, the Hipparcos satellite recently provided 
geometric parallax measurements for many nearby old stars with low metallicity, typical of glubular clus- 
ters, thus allowing for a new calibration of the ages of stars in globular clusters, leading to a downward 
revision to 10 — 13 Gyr 162) . Moreover, there were very few stars in the Hipparcos catalog with both 
small parallax erros and low metal abundance. Hence, an increase in the sample size could be critical in 
reducing the statatistical uncertainties for the calibration of the globular cluster ages. There are already 
proposed two new parallax satellites, NASA's Space Interferometry Mission (SIM) and ESAs mission, 
called GAIA, that will give 2 or 3 orders of magnitude more accurate parallaxes than Hipparcos, down 
to fainter magnitude limits, for several orders of magnitude more stars. Until larger samples are avail- 
able, however, distance errors are likely to be the largest source of systematic uncertainty to the globular 
cluster age \ 29\. 
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Fig. 22: The recent estimates of the age of the universe and that of the oldest objects in our galaxy. The last three points 
correspond to the combined analysis of 8 different measurements, for h — 0.64, 0.68 and 7.2, which indicates a relatively weak 
dependence on h. The age of the Sun is accurately known and is included for reference. Error bars indicate la limits. The 
averages of the ages of the Galactic Halo and Disk are shaded in gray. Note that there isn't a single age estimate more than 
2a away from the average. The result to > t ga i is logically inevitable, but the standard EdS model does not satisfy this unless 
h < 0.55. FromRef. 02). 



The supernovae groups can also determine the age of the universe from their high redshift observa- 
tions. The high confidence regions in the (Qm, ^a) plane are almost parallel to the contours of constant 
age. For any value of the Hubble constant less than Hq = 70 km/s/Mpc, the implied age of the universe 
is greater than 13 Gyr, allowing enough time for the oldest stars in globular clusters to evolve [62). In- 
tegrating over Qm and the best fit value of the age in Hubble-time units is H^q = 0.93 ± 0.06 or 
equivalently to = 14.1 ± 1.0 (0.65 h^ 1 ) Gyr, see Ref. [7 1. Furthermore, a combination of 8 independent 
recent measurements: CMB anisotropies, type la SNe, cluster mass-to-light ratios, cluster abundance 
evolution, cluster baryon fraction, deuterium-to-hidrogen ratios in quasar spectra, double-lobed radio 
sources and the Hubble constant, can be used to determine the present age of the universe [63|. The 



result is shown in Fig. |22l compared to other recent determinations. The best fit value for the age of the 
universe is, according to this analysis, to = 13.4 ± 1.6 Gyr, about a billion years younger than other 
recent estimates [63 1. 




Fig. 23: The anisotropies of the microwave background measured by the WMAP satellite with 10 arcminute resolution. It 
shows the intrinsic CMB anisotropies at the level of a few parts in 10 5 . The galactic foreground has been properly subtracted. 
The amount of information contained in this map is enough to determine most of the cosmological parameters to few percent 
accuracy. From Ref. |20|. 



3.4 Cosmic Microwave Background Anisotropies 

The cosmic microwave background has become in the last five years the Holy Grail of Cosmology, since 
precise observations of the temperature and polarization anisotropies allow in principle to determine the 
parameters of the Standard Model of Cosmology with very high accuracy. Recently, the WMAP satellite 
has provided with a very detailed map of the microwave anisotropies in the sky, see Fig. 23, and indeed 
has fulfilled our expectations, see Table 2. 

The physics of the CMB anisotropies is relatively simple [64|. The universe just before re- 
combination is a very tightly coupled fluid, due to the large electromagnetic Thomson cross section 
ot = 87ra 2 /3mg ~ 0.7 barn. Photons scatter off charged particles (protons and electrons), and carry 
energy, so they feel the gravitational potential associated with the perturbations imprinted in the metric 
during inflation. An overdensity of baryons (protons and neutrons) does not collapse under the effect of 
gravity until it enters the causal Hubble radius. The perturbation continues to grow until radiation pres- 
sure opposes gravity and sets up acoustic oscillations in the plasma, very similar to sound waves. Since 
overdensities of the same size will enter the Hubble radius at the same time, they will oscillate in phase. 
Moreover, since photons scatter off these baryons, the acoustic oscillations occur also in the photon field 
and induces a pattern of peaks in the temperature anisotropies in the sky, at different angular scales, see 
Fig. 24. There are three different effects that determine the temperature anisotropies we observe in the 
CMB. First, gravity: photons fall in and escape off gravitational potential wells, characterized by in the 
comoving gauge, and as a consequence their frequency is gravitationally blue- or red-shifted, 5v jv = 
If the gravitational potential is not constant, the photons will escape from a larger or smaller potential 
well than they fell in, so their frequency is also blue- or red-shifted, a phenomenon known as the Rees- 
Sciama effect. Second, pressure: photons scatter off baryons which fall into gravitational potential wells 
and the two competing forces create acoustic waves of compression and rarefaction. Finally, velocity: 
baryons accelerate as they fall into potential wells. They have minimum velocity at maximum compres- 
sion and rarefaction. That is, their velocity wave is exactly 90° off-phase with the acoustic waves. These 
waves induce a Doppler effect on the frequency of the photons. The temperature anisotropy induced by 
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Fig. 24: The Angular Power Spectrum of CMB temperature anisotropies, compared with the cross-correlation of temperature- 
polarization anisotropies. FromRef. |20|. 



these three effects is therefore given by l64l 
5T 
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Metric perturbations of different wavelengths enter the horizon at different times. The largest wave- 
lengths, of size comparable to our present horizon, are entering now. There are perturbations with wave- 
lengths comparable to the size of the horizon at the time of last scattering, of projected size about 1° 
in the sky today, which entered precisely at decoupling. And there are perturbations with wavelengths 
much smaller than the size of the horizon at last scattering, that entered much earlier than decoupling, all 
the way to the time of radiation-matter equality, which have gone through several acoustic oscillations 
before last scattering. All these perturbations of different wavelengths leave their imprint in the CMB 
anisotropies. 

The baryons at the time of decoupling do not feel the gravitational attraction of perturbations with 
wavelength greater than the size of the horizon at last scattering, because of causality. Perturbations with 
exactly that wavelength are undergoing their first contraction, or acoustic compression, at decoupling. 
Those perturbations induce a large peak in the temperature anisotropies power spectrum, see Fig. 24. 
Perturbations with wavelengths smaller than these will have gone, after they entered the Hubble scale, 
through a series of acoustic compressions and rarefactions, which can be seen as secondary peaks in 
the power spectrum. Since the surface of last scattering is not a sharp discontinuity, but a region of 
Az ~ 100, there will be scales for which photons, travelling from one energy concentration to another, 
will erase the perturbation on that scale, similarly to what neutrinos or HDM do for structure on small 
scales. That is the reason why we don't see all the acoustic oscillations with the same amplitude, but in 
fact they decay exponentialy towards smaller angular scales, an effect known as Silk damping, due to 
photon diffusion 1651 16411 . 



Table 2: The parameters of the standard cosmological model. The standard model of cosmology has about 20 different 
parameters, needed to describe the background space-time, the matter content and the spectrum of metric perturbations. We 
include here the present range of the most relevant parameters (with lu errors), as recently determined by WMAP, and the error 
with which the Planck satellite will be able to determine them in the near future. The rate of expansion is written in units of 
H = 100 h km/s/Mpc. 



physical quantity 


symbol 


WMAP 


Planck 


total density 




1.02 ±0.02 


0.7% 


baryonic matter 




0.044 ± 0.004 


0.6% 


cosmological constant 


n A 


0.73 ± 0.04 


0.5% 


cold dark matter 




0.23 ± 0.04 


0.6% 


hot dark matter 


n u h 2 


< 0.0076 (95% c.l.) 


1% 


sum of neutrino masses 


J2 m v (eV) 


< 0.23 (95% c.l.) 


1% 


CMB temperature 


To(K) 


2.725 ±0.002 


0.1% 


baryon to photon ratio 


i] x 10 iU 


6.1 ±0.3 


0.5% 


baryon to matter ratio 




0.17 ±0.01 


1% 


spatial curvature 


n K 


< 0.02 (95% c.l.) 


0.5% 


rate of expansion 


h 


0.71 ±0.03 


0.8% 


age of the universe 


t (Gyr) 


13.7 ±0.2 


0.1% 


age at decoupling 


idee (kyr) 


379 ±8 


0.5% 


age at reionization 


tr (Myr) 


180 ± 100 


5% 


spectral amplitude 


A 


0.833 ± 0.085 


0.1% 


spectral tilt 


n s 


0.98 ± 0.03 


0.2% 


spectral tilt variation 


dn s /d In k 


-0.031 ±0.017 


0.5% 


tensor-scalar ratio 


r 


< 0.71 (95% c.l.) 


5% 


reionization optical depth 


T 


0.17 ±0.04 


5% 


redshift of equality 




3233 ± 200 


5% 


redshift of decoupling 


Zdec 


1089 ± 1 


0.1% 


width of decoupling 


Azdec 


195 ±2 


1% 


redshift of reionization 


Z T 


20 ± 10 


2% 



From the observations of the CMB anisotropics it is possible to determine most of the parame- 
ters of the Standard Cosmological Model with few percent accuracy, see Table 2. However, there are 
many degeneracies between parameters and it is difficult to disentangle one from another. For instance, 
as mentioned above, the first peak in the photon distribution corresponds to overdensities that have un- 
dergone half an oscillation, that is, a compression, and appear at a scale associated with the size of the 
horizon at last scattering, about 1° projected in the sky today. Since photons scatter off baryons, they 
will also feel the acoustic wave and create a peak in the correlation function. The height of the peak 
is proportional to the amount of baryons: the larger the baryon content of the universe, the higher the 
peak. The position of the peak in the power spectrum depends on the geometrical size of the particle 
horizon at last scattering. Since photons travel along geodesies, the projected size of the causal horizon 
at decoupling depends on whether the universe is flat, open or closed. In a flat universe the geodesies 
are straight lines and, by looking at the angular scale of the first acoustic peak, we would be measuring 
the actual size of the horizon at last scattering. In an open universe, the geodesies are inward-curved tra- 
jectories, and therefore the projected size on the sky appears smaller. In this case, the first acoustic peak 
should occur at higher multipoles or smaller angular scales. On the other hand, for a closed universe, 
the first peak occurs at smaller multipoles or larger angular scales. The dependence of the position of 



the first acoustic peak on the spatial curvature can be approximately given by l pea k ~ 220 $7 , where 
Qo = Qm+^a = I — Qr- Present observations by WMAP and other experiments give Oo = 1.00±0.02 
at one standard deviation [20|. 
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Fig. 25: The (Qm, Ha) plane with the present data set of cosmological observations — the acceleration of the universe, the 
large scale structure and the CMB anisotropics — as well as the future determinations by SNAP and Planck of the fundamental 
parameters which define our Standard Model of Cosmology. 

The other acoustic peaks occur at harmonics of this, corresponding to smaller angular scales. Since 
the amplitude and position of the primary and secondary peaks are directly determined by the sound 
speed (and, hence, the equation of state) and by the geometry and expansion of the universe, they can 
be used as a powerful test of the density of baryons and dark matter, and other cosmological parameters. 
With the joined data from WMAP, VS A, CBI and ACBAR, we have rather good evidence of the existence 
of the second and third acoustic peaks, which confirms one of the most important predictions of inflation 
— the non-causal origin of the primordial spectrum of perturbations — , and rules out cosmological defects 
as the dominant source of structure in the universe [66]. Moreover, since the observations of CMB 
anisotropics now cover almost three orders of magnitude in the size of perturbations, we can determine 
the much better accuracy the value of the spectral tilt, n = 0.98 db 0.03, which is compatible with the 
approximate scale invariant spectrum needed for structure formation, and is a prediction of the simplest 
models of inflation. Soon after the release of data from WMAP, there was some expectation at the claim 
of a scale-dependent tilt. Nowadays, with better resolution in the linear matter power spectrum from 
SDSS [67], we can not conclude that the spectral tilt has any observable dependence on scale. 

The microwave background has become also a testing ground for theories of particle physics. In 



particular, it already gives stringent constraints on the mass of the neutrino, when analysed together with 
large scale structure observations. Assuming a flat ACDM model, the 2-sigma upper bounds on the sum 
of the masses of light neutrinos is J2 Wi/ < 1-0 eV for degenerate neutrinos (i.e. without a large hierachy 
between them) if we don't impose any priors, and it comes down to J2 m v < 0.6 eV if one imposes 
the bounds coming from the HST measurements of the rate of expansion and the supernova data on the 
present acceleration of the universe [68 1. The final bound on the neutrino density can be expressed as 
&u h 2 = J2 m v/93.2 eV < 0.01. In the future, both with Planck and with the Atacama Cosmology 
Telescope (ACT) we will be able to put constraints on the neutrino masses down to the 0.1 eV level. 

Moreover, the present data is good enough that we can start to put constraints on the models of in- 
flation that give rise to structure. In particular, multifield models of inflation predict a mixture of adiabatic 
and isocurvature perturbations, 9 and their signatures in the cosmic microwave background anisotropics 
and the matter power spectrum of large scale structure are specific and perfectly distinguishable. Nowa- 
days, thanks to precise CMB, LSS and SNIa data, one can put rather stringent limits on the relative 
fraction and correlation of the isocurvature modes to the dominant adiabatic perturbations [69|. 

We can summarize this Section by showing the region in parameter space where we stand nowa- 
days, thanks to the recent cosmological observations. We have plotted that region in Fig. 25. One could 
also superimpose the contour lines corresponding to equal io^o lines, as a cross check. It is extraordi- 
nary that only in the last few months we have been able to reduce the concordance region to where it 
stands today, where all the different observations seem to converge. There are still many uncertainties, 
mainly systematic; however, those are quickly decreasing and becoming predominantly statistical. In the 
near future, with precise observations of the anisotropics in the microwave background temperature and 
polarization anisotropics, thanks to Planck satellite, we will be able to reduce those uncertainties to the 
level of one percent. This is the reason why cosmologists are so excited and why it is claimed that we 
live in the Golden Age of Cosmology. 

4. THE INFLATIONARY PARADIGM 

The hot Big Bang theory is nowadays a very robust edifice, with many independent observational checks: 
the expansion of the universe; the abundance of light elements; the cosmic microwave background; a 
predicted age of the universe compatible with the age of the oldest objects in it, and the formation of 
structure via gravitational collapse of initially small inhomogeneities. Today, these observations are 
confirmed to within a few percent accuracy, and have helped establish the hot Big Bang as the preferred 
model of the universe. All the physics involved in the above observations is routinely tested in the 
laboratory (atomic and nuclear physics experiments) or in the solar system (general relativity). 

However, this theory leaves a range of crucial questions unanswered, most of which are initial 
conditions' problems. There is the reasonable assumption that these cosmological problems will be 
solved or explained by new physical principles at high energies, in the early universe. This assumption 
leads to the natural conclusion that accurate observations of the present state of the universe may shed 
light onto processes and physical laws at energies above those reachable by particle accelerators, present 
or future. We will see that this is a very optimistic approach indeed, and that there are many unresolved 
issues related to those problems. However, there might be in the near future reasons to be optimistic. 

4.1 Shortcomings of Big Bang Cosmology 

The Big Bang theory could not explain the origin of matter and structure in the universe; that is, the 
origin of the matter-antimatter asymmetry, without which the universe today would be filled by a uniform 
radiation continuosly expanding and cooling, with no traces of matter, and thus without the possibility 
to form gravitationally bound systems like galaxies, stars and planets that could sustain life. Moreover, 

'This mixture is generic, unless all the fields thermalize simultaneously at reheating, just after inflation, in which case the 
entropy perturbations that would give rise to the isocurvature modes disappear. 



the standard Big Bang theory assumes, but cannot explain, the origin of the extraordinary smoothness 
and flatness of the universe on the very large scales seen by the microwave background probes and the 
largest galaxy catalogs. It cannot explain the origin of the primordial density perturbations that gave rise 
to cosmic structures like galaxies, clusters and superclusters, via gravitational collapse; the quantity and 
nature of the dark matter that we believe holds the universe together; nor the origin of the Big Bang itself. 

A summary [ 10 1 of the problems that the Big Bang theory cannot explain is: 

• The global structure of the universe. 

- Why is the universe so close to spatial flatness? 

- Why is matter so homogeneously distributed on large scales? 

• The origin of structure in the universe. 

- How did the primordial spectrum of density perturbations originate? 

• The origin of matter and radiation. 

- Where does all the energy in the universe come from? 

- What is the nature of the dark matter in the universe? 

- How did the matter-antimatter asymmetry arise? 

• The initial singularity. 

- Did the universe have a beginning? 

- What is the global structure of the universe beyond our observable patch? 
Let me discuss one by one the different issues: 



4.11 The Flatness Problem 

The Big Bang theory assumes but cannot explain the extraordinary spatial flatness of our local patch of 
the universe. In the general FRW metric © the parameter K that characterizes spatial curvature is a free 
parameter. There is nothing in the theory that determines this parameter a priori. However, it is directly 
related, via the Friedmann equation jSj, to the dynamics, and thus the matter content, of the universe, 

K = — - — pa — Ha = — - — pa ( — — — ) . (119) 



3 r 3 r V Q 

We can therefore define a new variable, 



— 1 const. 
Q pa 2 



x = = , (120) 



whose time evolution is given by 



dx 

x' = — = (l + 3co)x, (121) 

where N = ln(a/aj) characterizes the number of e-folds of universe expansion (dN = Hdt) and 
where we have used Eq. Q for the time evolution of the total energy, pa 3 , which only depends on the 
barotropic ratio to. It is clear from Eq. (I121t that the phase-space diagram (x, x') presents an unstable 
critical (saddle) point at x = for lu > —1/3, i.e. for the radiation (u = 1/3) and matter (lo = 0) eras. 
A small perturbation from x = will drive the system towards x = ±oo. Since we know the universe 
went through both the radiation era (because of primordial nucleosynthesis) and the matter era (because 
of structure formation), tiny deviations from Q = 1 would have grown since then, such that today 

Qq — 1 / T; n \ 2 

£o = — p: = x in — (1 + z cq ) . (122) 

i Lq \ 1 eq / 

In order that today's value be in the range 0.1 < Qo < 1-2, or xq « 0(1), it is required that at, say, 
primordial nucleosynthesis (T NS ~ 10 6 T eq ) its value be 

Q(t NS ) = 1 ± 10- 15 , (123) 



which represents a tremendous finetuning. Perhaps the universe indeed started with such a peculiar 
initial condition, but it is epistemologically more satisfying if we give a fundamental dynamical reason 
for the universe to have started so close to spatial flatness. These arguments were first used by Robert 
Dicke in the 1960s, much before inflation. He argued that the most natural initial condition for the 
spatial curvature should have been the Planck scale curvature, = 6K/ lp , where the Planck length 
is lp = (hG/c 3 ) 1 / 2 = 1.62 x 10 -33 cm, that is, 60 orders of magnitude smaller than the present size 
of the universe, ao = 1.38 x 10 28 cm. A universe with this immense curvature would have collapsed 
within a Planck time, t P = (hG/c 5 ) 1 / 2 = 5.39 x 10" 44 s, again 60 orders of magnitude smaller than 
the present age of the universe, to = 4.1 x 10 17 s. Therefore, the flatness problem is also related to the 
Age Problem, why is it that the universe is so old and flat when, under ordinary circumstances (based on 
the fundamental scale of gravity) it should have lasted only a Planck time and reached a size of order the 
Planck length? As we will see, inflation gives a dynamical reason to such a peculiar initial condition. 

4.12 The Homogeneity Problem 

An expanding universe has particle horizons, that is, spatial regions beyond which causal communica- 
tion cannot occur. The horizon distance can be defined as the maximum distance that light could have 
travelled since the origin of the universe lfT5l . 

d H (t) = a(t) / —-H-^t), (124) 
Jo a(t') 

which is proportional to the Hubble scale. 10 For instance, at the beginning of nucleosynthesis the horizon 
distance is a few light-seconds, but grows linearly with time and by the end of nucleosynthesis it is a 
few light-minutes, i.e. a factor 100 larger, while the scale factor has increased only a factor of 10. The 
fact that the causal horizon increases faster, ~ t, than the scale factor, a ~ t 1//2 , implies that at any 
given time the universe contains regions within itself that, according to the Big Bang theory, were never 
in causal contact before. For instance, the number of causally disconnected regions at a given redshift z 
present in our causal volume today, g?h(£o) = a o> is 

N ^~{M^ i+zf2 - 0251 

which, for the time of decoupling, is of order AcD^dec) ~ 10 5 ^> 1. 

This phenomenon is particularly acute in the case of the observed microwave background. Infor- 
mation cannot travel faster than the speed of light, so the causal region at the time of photon decoupling 
could not be larger than du(tdec) ~ 3 x 10 5 light years across, or about 1° projected in the sky today. So 
why should regions that are separated by more than 1° in the sky today have exactly the same tempera- 
ture, to within 10 ppm, when the photons that come from those two distant regions could not have been 
in causal contact when they were emitted? This constitutes the so-called horizon problem, see Fig. |26l 
and was first discussed by Robert Dicke in the 1970s as a profound inconsistency of the Big Bang theory. 

4.2 Cosmological Inflation 

In the 1980s, a new paradigm, deeply rooted in fundamental physics, was put forward by Alan H. 
Guth [71 1, Andrei D. Linde [72| and others ll73l l74l l75l. to address these fundamental questions. Ac- 
cording to the inflationary paradigm, the early universe went through a period of exponential expansion, 
driven by the approximately constant energy density of a scalar field called the inflaton. In modern 
physics, elementary particles are represented by quantum fields, which resemble the familiar electric, 
magnetic and gravitational fields. A field is simply a function of space and time whose quantum oscil- 
lations are interpreted as particles. In our case, the inflaton field has, associated with it, a large potential 

10 For the radiation era, the horizon distance is equal to the Hubble scale. For the matter era it is twice the Hubble scale. 
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Fig. 26: Perhaps the most acute problem of the Big Bang theory is explaining the extraordinary homogeneity and isotropy of the 
microwave background, see Fig. llOl At the time of decoupling, the volume that gave rise to our present universe contained many 
causally disconnected regions (top figure). Today we observe a blackbody spectrum of photons coming from those regions and 
they appear to have the same temperature, T\ = T%, to one part in 10 5 . Why is the universe so homogeneous? This constitutes 
the so-called horizon problem, which is spectacularly solved by inflation. From Ref. |70|. 

energy density, which drives the exponential expansion during inflation, see Fig.|^] We know from gen- 
eral relativity that the density of matter determines the expansion of the universe, but a constant energy 
density acts in a very peculiar way: as a repulsive force that makes any two points in space separate at 
exponentially large speeds. (This does not violate the laws of causality because there is no information 
carried along in the expansion, it is simply the stretching of space-time.) 

This superluminal expansion is capable of explaining the large scale homogeneity of our observ- 
able universe and, in particular, why the microwave background looks so isotropic: regions separated 
today by more than 1° in the sky were, in fact, in causal contact before inflation, but were stretched to 
cosmological distances by the expansion. Any inhomogeneities present before the tremendous expansion 
would be washed out. This explains why photons from supposedly causally disconneted regions have 
actually the same spectral distribution with the same temperature, see Fig. |26] 

Moreover, in the usual Big Bang scenario a flat universe, one in which the gravitational attraction 
of matter is exactly balanced by the cosmic expansion, is unstable under perturbations: a small deviation 
from flatness is amplified and soon produces either an empty universe or a collapsed one. As we dis- 
cussed above, for the universe to be nearly flat today, it must have been extremely flat at nucleosynthesis, 
deviations not exceeding more than one part in 10 15 . This extreme fine tuning of initial conditions was 
also solved by the inflationary paradigm, see Fig. |28] Thus inflation is an extremely elegant hypothesis 
that explains how a region much, much greater that our own observable universe could have become 
smooth and flat without recourse to ad hoc initial conditions. Furthermore, inflation dilutes away any 
"unwanted" relic species that could have remained from early universe phase transitions, like monopoles, 
cosmic strings, etc., which are predicted in grand unified theories and whose energy density could be so 
large that the universe would have become unstable, and collapsed, long ago. These relics are diluted by 
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Fig. 27: The inflaton field can be represented as a ball rolling down a hill. During inflation, the energy density is approximately 
constant, driving the tremendous expansion of the universe. When the ball starts to oscillate around the bottom of the hill, 
inflation ends and the inflaton energy decays into particles. In certain cases, the coherent oscillations of the inflaton could 
generate a resonant production of particles which soon thermalize, reheating the universe. From Ref. |70|. 



the superluminal expansion, which leaves at most one of these particles per causal horizon, making them 
harmless to the subsequent evolution of the universe. 

The only thing we know about this peculiar scalar field, the inflaton, is that it has a mass and 
a self-interaction potential V(4>) but we ignore everything else, even the scale at which its dynamics 
determines the superluminal expansion. In particular, we still do not know the nature of the inflaton field 
itself, is it some new fundamental scalar field in the electroweak symmetry breaking sector, or is it just 
some effective description of a more fundamental high energy interaction? Hopefully, in the near future, 
experiments in particle physics might give us a clue to its nature. Inflation had its original inspiration in 
the Higgs field, the scalar field supposed to be responsible for the masses of elementary particles (quarks 
and leptons) and the breaking of the electroweak symmetry. Such a field has not been found yet, and its 
discovery at the future particle colliders would help understand one of the truly fundamental problems in 
physics, the origin of masses. If the experiments discover something completely new and unexpected, it 
would automatically affect the idea of inflation at a fundamental level. 

4.21 Homogeneous scalar field dynamics 

In this subsection I will describe the theoretical basis for the phenomenon of inflation. Consider a scalar 
field <p, a singlet under any given interaction, with an effective potential V(<p). The Lagrangian for such 
a field in a curved background is 

Anf ^sTd^t-VM, (126) 

whose evolution equation in a Friedmann-Robertson-Walker metric © and for a homogeneous field cj)(t) 
is given by 

4> + 3H<P + V'{(f>) = , (127) 
where H is the rate of expansion, together with the Einstein equations, 

H 2 = y(^ 2 + V^)), (128) 



Fig. 28: The exponential expansion during inflation made the radius of curvature of the universe so large that our observable 
patch of the universe today appears essentialy flat, analogous (in three dimensions) to how the surface of a balloon appears 
flatter and flatter as we inflate it to enormous sizes. This is a crucial prediction of cosmological inflation that will be tested to 
extraordinary accuracy in the next few years. From Ref. 1741 1701 . 



H = -y^ 2 , (129) 

where k 2 = 8ttG. The dynamics of inflation can be described as a perfect fluid © with a time dependent 
pressure and energy density given by 

9 = \4> 2 + V{^), (130) 

P = ^ 2 -V(<p). (131) 

The field evolution equation dl27t can then be written as the energy conservation equation, 

p + ZH{p + p) = . (132) 

If the potential energy density of the scalar field dominates the kinetic energy, V{(f>) S> <fi 2 , then we see 
that 

p ~ — p p ~ const. =^ H(4>) ~ const. , (133) 

which leads to the solution 

a(t) ~ exp(Ht) — >0 accelerated expansion . (134) 

a 

Using the definition of the number of e-folds, N = ln(a/oj), we see that the scale factor grows expo- 
nentially, a(N) = Ojexp(iV). This solution of the Einstein equations solves immediately the flatness 



problem. Recall that the problem with the radiation and matter eras is that Q, = 1 (x = 0) is an un- 
stable critical point in phase-space. However, during inflation, with p ~ — p =4> u ~ —1, we have 
that 1 + 3lo > and therefore x = is a stable attractor of the equations of motion, see Eq. (I12H . 
As a consequence, what seemed an ad hoc initial condition, becomes a natural prediction of inflation. 
Suppose that during inflation the scale factor increased N e-folds, then 

X0 = Xine -2^(^l) 2 ( i + Zcq) ~ e - 2N W™<l N> 65, (135) 

-*eq 

where we have assumed that inflation ended at the scale V en d, and the transfer of the inflaton energy 
density to thermal radiation at reheating occurred almost instantaneously 11 at the temperature ~ 

Knd ~ 1Ql5 GeV - N ° te that 

we can now have initial conditions with a large uncertainty, x- m ~ 1, and 
still have today xq ~ 1, thanks to the inflationary attractor towards Q, = 1. This can be understood very 
easily by realizing that the three curvature evolves during inflation as 

® R = ^£ = m R . n e -2N _^ for iV » 1 . (136) 

Therefore, if cosmological inflation lasted over 65 e-folds, as most models predict, then today the uni- 
verse (or at least our local patch) should be exactly flat, see Fig. |28j a prediction that can be tested with 
great accuracy in the near future and for which already seems to be some evidence from observations of 
the microwave background [ 87 1. 

Furthermore, inflation also solves the homogeneity problem in a spectacular way. First of all, due 
to the superluminal expansion, any inhomogeneity existing prior to inflation will be washed out, 



(~Ir) $fc K e ~ m ~~ * ' for N > 1 ' 



(137) 



Moreover, since the scale factor grows exponentially, while the horizon distance remains essentially 
constant, djj{t) ~ H^ 1 = const., any scale within the horizon during inflation will be stretched by the 
superluminal expansion to enormous distances, in such a way that at photon decoupling all the causally 
disconnected regions that encompass our present horizon actually come from a single region during 
inflation, about 65 e-folds before the end. This is the reason why two points separated more than 1° 
in the sky have the same backbody temperature, as observed by the COBE satellite: they were actually 
in causal contact during inflation. There is at present no other proposal known that could solve the 
homogeneity problem without invoquing an acausal mechanism like inflation. 

Finally, any relic particle species (relativistic or not) existing prior to inflation will be diluted by 
the expansion, 

pu oc a" 3 ~ e~ m — > 0, foriV>l, (138) 
Pr oc a" 4 ~ e~ 4N — ► 0, foriV>l. (139) 

Note that the vacuum energy density p v remains constant under the expansion, and therefore, very soon 
it is the only energy density remaining to drive the expansion of the universe. 



4.22 The slow-roll approximation 

In order to simplify the evolution equations during inflation, we will consider the slow-roll approximation 
(SRA). Suppose that, during inflation, the scalar field evolves very slowly down its effective potential, 

"There could be a small delay in thermalization, due to the intrinsic inefficiency of reheating, but this does not change 
significantly the required number of e-folds. 



then we can define the slow -roll parameters [76|, 



H k 2 2 

^ = y^« 1 ' (140) 



Hcj) 



« 1, (141) 



It is easy to see that the condition 



£ = — U-<5 2 < 1. (142) 

e < 1 ^ - > (143) 

a 

characterizes inflation: it is all you need for superluminal expansion, i.e. for the horizon distance to grow 
more slowly than the scale factor, in order to solve the homogeneity problem, as well as for the spatial 
curvature to decay faster than usual, in order to solve the flatness problem. 

The number of e-folds during inflation can be written with the help of Eq. d!40t as 

iV = ln^ = ( tB Hdt= (144) 

which is an exact expression in terms of e{4>). 

In the limit given by Eqs. (I140I) . the evolution equations (I127t and dl28t become 

H 2 (l- € -) * H 2 = ^-V(cP), (145) 

3Hj>(l--) ~ 3H<j> = - V {<!>). (146) 

Note that this corresponds to a reduction of the dimensionality of phase-space from two to one dimen- 
sions, H((f>, 4>) — > H(4>)- In fact, it is possible to prove a theorem, for single-field inflation, which 
states that the slow-roll approximation is an attractor of the equations of motion, and thus we can al- 
ways evaluate the inflationary trajectory in phase-space within the SRA, therefore reducing the number 
of initial conditions to just one, the initial value of the scalar field. If H{(j>) only depends on 0, then 
H'(<p) = —K 2 cf)/2 and we can rewrite the slow-roll parameters dl40t as 



2 -, ,fTlt ±\ \ 2 

= e v « 1 , (147) 



2_(HWY „ 1 (V'(<t>)\ 
k 2 \ H{<t>) ) ~ 2k 2 \ V(</>) J 

2 h"(4>) i v"(4>) i fv'ty) 



k 2 H(4>) k 2 V{<f>) 2k 2 V V(</>) 
f 4 H'^H'"^) 1 ng^j 3 V"(d>) (V'{<P)\ 

1 H 2 (<p) k 4 V 2 (<p) 2k 4 V(4>) \V(<j>) J 



W-ev <C 1, (148) 

+ Mm)' Ber -*" cr+ *** 1 - <149) 

These expressions define the new slow -roll parameters ey, r/y and £y. The number of e-folds can also 
be rewritten in this approximation as 



= (150) 



a very useful expression for evaluating N for a given effective scalar potential V (<f>). 



4.3 The origin of density perturbations 

If cosmological inflation made the universe so extremely flat and homogeneous, where did the galaxies 
and clusters of galaxies come from? One of the most astonishing predictions of inflation, one that was not 
even expected, is that quantum fluctuations of the inflaton field are stretched by the exponential expansion 
and generate large-scale perturbations in the metric. Inflaton fluctuations are small wave packets of 
energy that, according to general relativity, modify the space-time fabric, creating a whole spectrum of 
curvature perturbations. The use of the word spectrum here is closely related to the case of light waves 
propagating in a medium: a spectrum characterizes the amplitude of each given wavelength. In the 
case of inflation, the inflaton fluctuations induce waves in the space-time metric that can be decomposed 
into different wavelengths, all with approximately the same amplitude, that is, corresponding to a scale- 
invariant spectrum. These patterns of perturbations in the metric are like fingerprints that unequivocally 
characterize a period of inflation. When matter fell in the troughs of these waves, it created density 
perturbations that collapsed gravitationally to form galaxies, clusters and superclusters of galaxies, with 
a spectrum that is also scale invariant. Such a type of spectrum was proposed in the early 1970s (before 
inflation) by Harrison and Zel'dovich [26 ], to explain the distribution of galaxies and clusters of galaxies 
on very large scales in our observable universe. Perhaps the most interesting aspect of structure formation 
is the possibility that the detailed knowledge of what seeded galaxies and clusters of galaxies will allow 
us to test the idea of inflation. 



4.31 Reparametrization invariant perturbation theory 

Until now we have considered only the unperturbed FRW metric described by a scale factor a(t) and a 
homogeneous scalar field <p(t), 

ds 2 = a 2 {rj)[-dr) 2 + j ij dx i dx j ], (151) 
<f> = <Kv), (152) 

where rj = f dt/a(t) is the conformal time, under which the background equations of motion can be 
written as 



n 2 = y (~0 ,2 + «M0)) , (153) 



H'-H 2 = -y0' 2 , (154) 

(/)" + 2Hcp' + a 2 V'{4)) = 0, (155) 

where H = aH and cf>' = acf>. 

During inflation, the quantum fluctuations of the scalar field will induce metric perturbations which 
will backreact on the scalar field. Let us consider, in linear perturbation theory, the most general line ele- 
ment with both scalar and tensor metric perturbations [77 ], 12 together with the scalar field perturbations 

ds 2 = a 2 (r]) [- (1 + 2A)dr] 2 + IB^dr) + {(1 + 211)^ + 2E {ij + 2h ij }dx i dx j ^ , (156) 
= ( t>( v )+S(l>(r ) ,x i ). (157) 

The indices {i,j} label the three-dimensional spatial coordinates with metric 7^, and the \i denotes 
covariant derivative with respect to that metric. The gauge invariant tensor perturbation h%j corresponds 
to a transverse traceless gravitational wave, V l hij = h\ = 0. The four scalar perturbations (A, B, TZ, E) 
are gauge dependent functions of (77, x l ). Under a general coordinate (gauge) transformation 17711781 

r? = 77 + £° (??,<), (158) 
5 i = x i + 7 %(^s i ), (159) 



''Note that inflation cannot generate, to linear order, a vector perturbation. 



with arbitrary functions (£ ,£), the scalar and tensor perturbations transform, to linear order, as 

A = A - f - , B = B + f , (160) 
K = K-HZ°, E = E-i, (161) 

hi] = hij , (162) 

where a prime denotes derivative with respect to conformal time. It is possible to construct, however, 
two gauge-invariant gravitational potentials [77 78 1, 



$ = A+ (B - E')' + H(B - E') , (163) 
= K + H(B-E'), (164) 



which are related through the perturbed Einstein equations, 



$ = (165) 

— = — Sp, (166) 

a z 2 

where 5 pis the gauge-invariant density perturbation, and the latter expression is nothing but the Poisson 
equation for the gravitational potential, written in relativistic form. 

During inflation, the energy density is given in terms of a scalar field, and thus the gauge-invariant 
equations for the perturbations on comoving hypersurfaces (constant energy density hypersurfaces) are 

K 2 

<5>" + m& + (H' + 2H 2 )§ = —[(f>'54>'-a 2 V'{4>)54>\, (167) 



K 2 

-V 2 $ + 3W$' + (W + 2H 2 )$> = - — \<t>'5<j)' + a 2 V'(<f))d(f>\, (168) 

k 2 

$' + H<5> = —<f>'5<t>, (169) 
5c/)" + 2H5(f)' -V 2 5<p = 4<j>'& - 2a 2 V {(f))<5> - a 2 V" {<]))5<]) . (170) 

This system of equations seem too difficult to solve at first sight. However, there is a gauge 
invariant combination of variables that allows one to find exact solutions. Let us define 1751 

u = a5(f) + z^, (171) 

z = a%. (172) 
Ft 

Under this redefinition, the above equations simplify enormously to just three independent equations, 



ji 



u" - V 2 u- — u = 0, (173) 
V 2 $ = ^- (zu - z'u) , (174) 



z 



2 a 2 ' 



/a 2 §\i k 

(-Jf)=Y zu - (175) 



From Equation d 173l > we can find a solution u(z), which substituted into dl75t can be integrated to give 
$>(z), and together with u(z) allow us to obtain 5<j){z). 



4.32 Quantum Field Theory in curved space-time 

Until now we have treated the perturbations as classical, but we should in fact consider the perturbations 
$ and 5(p as quantum fields. Note that the perturbed action for the scalar mode u can be written as 



1 



5S = - I d A xdr) {u'f - {Vuf + —u 



(176) 



In order to quantize the field u in the curved background defined by the metric (1151 b . we can write the 
operator 

d 3 k 



u(r),x) 



Wfc(j?)a k e l ' x + u* k {r{) a\,e 



t „-!k x 



(177) 



(2vr) 3 / 2 

where the creation and annihilation operators satisfy the commutation relation of bosonic fields, and the 
scalar field's Fock space is defined through the vacuum condition, 

[«k,4'l = 5 3 (k-k'), (178) 
o k |0> = 0. (179) 

Note that we are not assuming that the inflaton is a fundamental scalar field, but that is can be written as 
a quantum field with its commutation relations (as much as a pion can be described as a quantum field). 

The equations of motion for each mode u^{rj) are decoupled in linear perturbation theory, 



€ + ( fc2 - — ) u k = • 



(180) 



The ratio z" / z acts like a time-dependent potential for this Schrodinger like equation. In order to find 
exact solutions to the mode equation, we will use the slow-roll parameters dl40t . see Ref. ll76l 



K 2 Z 2 
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f 



1 + e 



Hz ' 



2 - e - 35 + 5 2 



■H 2 4>> 



(181) 
(182) 
(183) 



In terms of these parameters, the conformal time and the effective potential for the u k mode can be 
written as 



-1 



V 



z 
z 



n 



+ 



eda 



{l + e-5){2-5)+H- l {e' -8') . 
Note that the slow -roll parameters, (11811) and (I182I) . can be taken as constant, 13 to order e 2 , 

e' = 2H (e 2 - e6) = 0{e 2 ) , 
5> = H (eS -C)= 0(e 2 ) . 

In that case, for constant slow-roll parameters, we can write 

-1 1 



z rf V 4/ 



n l-e 
where v 



1+e 



6 1 
l- e - + 2 



(184) 
(185) 

(186) 

(187) 
(188) 



13 For instance, there are models of inflation, like power-law inflation, a(t) ~ t p , where e = 5 — 1/p < 1, that give constant 
slow-roll parameters. 



We are now going to search for approximate solutions of the mode equation d 1 801) . where the 
effective potential (11851) is of order z" / z ~ 2TL 2 in the slow-roll approximation. In quasi-de Sitter there 
is a characteristic scale given by the (event) horizon size or Hubble scale during inflation, H . There 
will be modes u k with physical wavelengths much smaller than this scale, k/a S> H, that are well 
within the de Sitter horizon and therefore do not feel the curvature of space-time. On the other hand, 
there will be modes with physical wavelengths much greater than the Hubble scale, k/a <C H. In these 
two asymptotic regimes, the solutions can be written as 

u k = -= e- ikri k^>aH , (189) 



f 2k 

u k = dz k<.aH . (190) 

In the limit k S> aH the modes behave like ordinary quantum modes in Minkowsky space-time, ap- 
propriately normalized, while in the opposite limit, u/z becomes constant on superhorizon scales. For 
approximately constant slow-roll parameters one can find exact solutions to d 1 801) . with the effective 
potential given by (11881) . that interpolate between the two asymptotic solutions, 

uifcfa) = ^ j {v+m (~V) 1/2 HP (-k v ) , (191) 

where Hy ] (z) is the Hankel function of the first kind 1791 . and v is given by (11881) in terms of the 
slow-roll parameters. In the limit k-q — ► 0, the solution becomes 

2"- a T(u) . , ,i u Civ) ( k \\-v 

^vmM^ = ^eU • (192) 



3 



r(l) 



C{v) = 2 U ~2 (1 - e)"-2 ~ i for e,S<l. (193) 



We can now compute $ and 5(p from the super-Hubble-scale mode solution ( I190I ). for k <C aH. 
Substituting into Eq. ill 51 . we find 

$ = Ci (l - ^ J a 2 dr]) +C 2 ^, (194) 
/ a 2 dr? - -| . (195) 



a* 

The term proportional to Ci corresponds to the growing solution, while that proportional to C 2 corre- 
sponds to the decaying solution, which can soon be ignored. These quantities are gauge invariant but 
evolve with time outside the horizon, during inflation, and before entering again the horizon during the 
radiation or matter eras. We would like to write an expression for a gauge invariant quantity that is 
also constant for superhorizon modes. Fortunately, in the case of adiabatic perturbations, there is such a 
quantity: 

1 u 
C = <D + — + = (196) 
en z 

which is constant, see Eq. (11901) . for k <C aH. In fact, this quantity £ is identical, for superhorizon 
modes, to the gauge invariant curvature metric perturbation 1Z C on comoving (constant energy density) 
hypersurfaces, see Ref. [77, 80], 

C = ^ C + ^V 2 $. (197) 

Using Eq. d!74t we can write the evolution equation for ( = ~ as (' = i V 2 <&, which confirms that 
£ is constant for (adiabatic 14 ) superhorizon modes, k <C aH. Therefore, we can evaluate the Newtonian 

14 This conservation fails for entropy or isocurvature perturbations, see Ref. 1801 . 



potential & k when the perturbation reenters the horizon during radiation/matter eras in terms of the 
curvature perturbation 1Z k when it left the Hubble scale during inflation, 

n-i r \ o _i_ q, , f | IZu radiation era , 

« 2 dv )K k = l±^K k = (198) 
v a J ' b + du [ matter era. 

Let us now compute the tensor or gravitational wave metric perturbations generated during infla- 
tion. The perturbed action for the tensor mode can be written as 



1 r 43 - a 



2 



6S= 2j d x dv 2^ ~ 1 ' (199) 

with the tensor field h{j considered as a quantum field, 

f d 3 k r l- i 

hij(ri,x) = J ^ 3/2 [M^) M k > A ) «k,A e lkoc + /i.c.J , (200) 

where &ij(k, A) are the two polarization tensors, satisfying symmetric, transverse and traceless conditions 

&ij — ^ji 7 ^ ^ij — j &ii — , (201) 

ey(-k, A) = e*-(k, A) , £ e*(k, X)e ij (k, A) = 4 , (202) 

A 

while the creation and annihilation operators satisfy the usual commutation relation of bosonic fields, 
Eq. (11781 . We can now redefine our gauge invariant tensor amplitude as 

v k (v) = -%-h k (r t ), (203) 

which satisfies the following evolution equation, decoupled for each mode v k (rj) in linear perturbation 
theory, 

vl+ (k 2 -— )v k = 0. (204) 



a 

The ratio a" /a acts like a time-dependent potential for this Schrodinger like equation, analogous to the 
term z" / z for the scalar metric perturbation. For constant slow-roll parameters, the potential becomes 

-£)_£(„>-!), 

f=ih + l (206) 

We can solve equation (12041) in the two asymptotic regimes, 

v k = -== e~ ikr] k^>aH, (207) 
y/2k 



v k = Ca k<^aH . (208) 

In the limit k 3> aH the modes behave like ordinary quantum modes in Minkowsky space-time, ap- 
propriately normalized, while in the opposite limit, the metric perturbation h k becomes constant on 
superhorizon scales. For constant slow-roll parameters one can find exact solutions to (1204 1 . with effec- 
tive potential given by d205t . that interpolate between the two asymptotic solutions. These are identical 
to Eq. d!91t except for the substitution v — ► fi. In the limit krj — ► 0, the solution becomes 



Since the mode h k becomes constant on superhorizon scales, we can evaluate the tensor metric pertur- 
bation when it reentered during the radiation or matter era directly in terms of its value during inflation. 



4.33 Power spectrum of scalar and tensor metric perturbations 

Not only do we expect to measure the amplitude of the metric perturbations generated during inflation 
and responsible for the anisotropics in the CMB and density fluctuations in LSS, but we should also be 
able to measure its power spectrum, or two-point correlation function in Fourier space. Let us consider 
first the scalar metric perturbations TZf,, which enter the horizon at a = k/H. Its correlator is given 
by UB 

m * k K k ,\0) = \^5 3 (k- k') ee (2tt) 3 5 3 (k - k') , (210) 

<d fu\ k3 l Mfc l 2 ( H \ 2 ( k \ 3 ~ 2u - a2 ( k V s - 1 n\\\ 

V ^ k) = 2^— = Te^) Krf) = M^J ' (2H) 

where we have used TZk = (k = an d Eq. (I192I) . This last equation determines the power spectrum in 
terms of its amplitude at horizon-crossing, As, and a tilt, 

dlnP n (k) „ n f6-2e^ 



n s ~l = - 11 = [(1 - e) - (e - 5) - l] = 2^j— ^) ^ 2rj v - Qe v , (213) 



1 - ^PF = 3 " 2u = 2(°-^) * 2 Vv - 6e v , (212) 
a In A; M-e' 

see Eqs. d 1 47 1 ) . d!48t . Note from this equation that it is possible, in principle, to obtain from inflation a 
scalar tilt which is either positive (n > 1) or negative (n < 1). Furthermore, depending on the particular 
inflationary model [81 1, we can have significant departures from scale invariance. 

Note that at horizon entry ki] = — 1, and thus we can alternatively evaluate the tilt as 

'S-2e 

- - ~->>iH i I — ( ) — (( - ,)) - I - •>[ 

din 7] 

and the running of the tilt 

dhk = ~ dh^ = V 2 ^ + 8e2 " We V ~ 2 ^ v + 24e ^ " 1671vev ' (214) 

where we have used Eqs. d!86t . 

Let us consider now the tensor (gravitational wave) metric perturbation, which enter the horizon 

at a = k/H, 

E<°I*M VaIO) = 4 ^ |^| 2 5 3 (k - k') ee ^ (2vr) 3 <5 3 (k - k') , (215) 

A 



where we have used Eqs. (1203 1 and (1209b . Therefore, the power spectrum can be approximated by a 
power-law expression, with amplitude At and tilt 

^ ^ dJnP^ = g -2e_ ^ ^ 

dink 1 — e 

which is always negative. In the slow-roll approximation, e -C 1, the tensor power spectrum is scale 
invariant. 

Alternatively, we can evaluate the tensor tilt by 



dhxV, 

tit ee — 

and its running by 



nT = _ ^±L£ = _2r]H [(1 - e) - ll = — ?i ~ -2e v , (218) 
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where we have used Eqs. 
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4.4 The anisotropics of the microwave background 

The metric fluctuations generated during inflation are not only responsible for the density perturbations 
that gave rise to galaxies via gravitational collapse, but one should also expect to see such ripples in 
the metric as temperature anisotropies in the cosmic microwave background, that is, minute deviations 
in the temperature of the blackbody spectrum when we look at different directions in the sky. Such 
anisotropies had been looked for ever since Penzias and Wilson's discovery of the CMB, but had eluded 
all detection, until COBE satellite discovered them in 1992, see Fig. ^0] The reason why they took so 
long to be discovered was that they appear as perturbations in temperature of only one part in 10 5 . Soon 
after COBE, other groups quickly confirmed the detection of temperature anisotropies at around 30 pK, 
at higher multipole numbers or smaller angular scales. 



4.41 The Sachs-Wolfe effect 

The anisotropies corresponding to large angular scales are only generated via gravitational red-shift 
and density perturbations through the Einstein equations, Sp/p = —2$ for adiabatic perturbations; we 
can ignore the Doppler contribution, since the perturbation is non-causal. In that case, the temperature 
anisotropy in the sky today is given by [ 82 1 
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where t]q is the coordinate distance to the last scattering surface, i.e. the present conformal time, while 
Vls — determines that comoving hypersurface. The above expression is known as the Sachs-Wolfe 
effect [82], and contains two parts, the intrinsic and the Integrated Sachs-Wolfe (ISW) effect, due to 
integration along the line of sight of time variations in the gravitational potential. 

In linear perturbation theory, the scalar metric perturbations can be separated into $(jj,x) = 
$(?7) Q(x), where Q( x ) are the scalar harmonics, eigenfunctions of the Laplacian in three dimensions, 
V 2 Qfe/ m (r, 6, (j)) = —k 2 Qkim{f, 0, (ft). These functions have the general form [83] 

Qki m (r, 9, 4>) = Il kl (r) Y lm (6, 0) , (221) 

where Yi m (8, <j)) are the usual spherical harmonics [79 1. 

In order to compute the temperature anisotropy associated with the Sachs-Wolfe effect, we have 
to know the evolution of the metric perturbation during the matter era, 

+ m $' + a 2 A $ - 2K 3> = . (222) 

In the case of a flat universe without cosmological constant, the Newtonian potential remains constant 
during the matter era and only the intrinsic SW effect contributes to 5T/T. In case of a non- vanishing A, 
since its contribution is negligible in the past, most of the photon's trajectory towards us is unperturbed, 
and the only difference with respect to the A = case is an overall factor [ 86 1. We will consider here the 
approximation $ = constant during the matter era and ignore that factor, see Ref. [84|. 

In a flat universe, the radial part of the eigenfunctions (122 It can be written as [ 83 1 

U kl (r) = ^k Jl (kr), (223) 

V 7T 

where ji{z) are the spherical Bessel functions [ 79 1 . The growing mode solution of the metric perturbation 
that left the Hubble scale during inflation contributes to the temperature anisotropies on large scales (I220t 
as 

ST 1 1 oa l 

— (0,0) = -$(77ls)Q = - ^ Q(r? , = EE a lm Y lm (e,cf>), (224) 
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where we have used the fact that at reentry (at the surface of last scattering) the gauge invariant Newtonian 
potential <3? is related to the curvature perturbation 1Z at Hubble-crossing during inflation, see Eq. ( I198I ): 
and we have expanded 5T/T in spherical harmonics. 

We can now compute the two-point correlation function or angular power spectrum, C{6), of the 
CMB anisotropics on large scales, defined as an expansion in multipole number, 

C(9) = (— (n) — (n')) = 5- £( 2 * + l)QPi(cos9) , (225) 
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where Pi(z) are the Legendre polynomials [79|, and we have averaged over different universe realiza- 
tions. Since the coefficients ai m are isotropic (to first order), we can compute the C\ = (|a; m | 2 ) as 

<S) _ 4vr [°° dk 2 
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where we have used Eqs. (12241) and d210t . In the case of scalar metric perturbation produced during 
inflation, the scalar power spectrum at reentry is given by Vn(k) = A 2 s (krio) n ~ l , in the power-law 
approximation, see Eq. (121 II) . In that case, one can integrate d226t to give 

(S ) 2n 9 rrfirri-Mirri + M] 

C} b} = — Al o i -i , (227) 

in + 1) d s) Al 

= — S = constant , for n = 1 . (228) 
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This last expression corresponds to what is known as the Sachs-Wolfe plateau, and is the reason why the 
coefficients Ci are always plotted multiplied by 1(1 + 1), see Fig. 13.41 

Tensor metric perturbations also contribute with an approximately constant angular power spec- 
trum, 1(1 + l)C[. The Sachs-Wolfe effect for a gauge invariant tensor perturbation is given by [82| 
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/ drh'( m -r)Q rr (r : 9,(t>), (229) 



where Q rr is the rr-component of the tensor harmonic along the line of sight [83 1. The tensor perturba- 
tion h during the matter era satisfies the following evolution equation 

h'{ + m h' k + (k 2 + 2K)h k = 0, (230) 

which depends on the wavenumber k, contrary to what happens with the scalar modes, see Eq. (1222b . For 
a flat (K = 0) universe, the solution to this equation is hk(fj) = h Gk(rj), where h is the constant tensor 
metric perturbation at horizon crossing and G k (rf) = 3ji(krj)/kr], normalized so that G k (0) = 1 at the 
surface of last scattering. The radial part of the tensor harmonic Q rr in a flat universe can be written 
as l83l 

"(Z -!)/(/ + !)(/ + 2)1 1/2 
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The tensor angular power spectrum can finally be expressed as 
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C (T) = JL(l- m + + 2) J" - V g (k) II , (232) 
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where x = krj, and V g (k) is the primordial tensor spectrum d2 16l >- For a scale invariant spectrum, 
ut = 0, we can integrate d232t to give [ 85 1 
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with Bi = (1.1184, 0.8789, . . . , 1.00) for I = 2, 3, ... , 30. Therefore, 1(1 + 1) C\ ' also becomes 
constant for large I. Beyond I ~ 30, the Sachs-Wolfe expression is not a good approximation and the 
tensor angular power spectrum decays very quickly at large I, see Fig.l3T1 

4.42 The consistency relation 

In spite of the success of inflation in predicting a homogeneous and isotropic background on which to 
imprint a scale-invariant spectrum of inhomogeneities, it is difficult to test the idea of inflation. A CMB 
cosmologist before the 1980s would have argued that ad hoc initial conditions could have been at the 
origin of the homogeneity and flatness of the universe on large scales, while a LSS cosmologist would 
have agreed with Harrison and Zel'dovich that the most natural spectrum needed to explain the formation 
of structure was a scale-invariant spectrum. The surprise was that inflation incorporated an understanding 
of both the globally homogeneous and spatially flat background, and the approximately scale-invariant 
spectrum of perturbations in the same formalism. But that could have been just a coincidence. 

What is unique to inflation is the fact that inflation determines not just one but two primordial 
spectra, corresponding to the scalar (density) and tensor (gravitational waves) metric perturbations, from 
a single continuous function, the inflaton potential V(4>). In the slow-roll approximation, one determines, 
from V(4>), two continuous functions, Vyi(k) and V g (k), that in the power-law approximation reduces 
to two amplitudes, As and At, and two tilts, n and tit- It is clear that there must be a relation between 
the four parameters. Indeed, one can see from Eqs. (I234t and (122 8 1 that the ratio of the tensor to scalar 
contribution to the angular power spectrum is proportional to the tensor tilt |76|, 



This is a unique prediction of inflation, which could not have been postulated a priori by any cosmol- 
ogist. If we finally observe a tensor spectrum of anisotropics in the CMB, or a stochastic gravitational 
wave background in laser interferometers like LIGO or LISA, with sufficient accuracy to determine their 
spectral tilt, one might have some chance to test the idea of inflation, via the consistency relation (I235t . 
For the moment, observations of the microwave background anisotropics suggest that the Sachs-Wolfe 
plateau exists, see Fig. 13.41 but it is still premature to determine the tensor contribution. Perhaps in the 
near future, from the analysis of polarization as well as temperature anisotropies, with the CMB satellites 
MAP and Planck, we might have a chance of determining the validity of the consistency relation. 

Assuming that the scalar contribution dominates over the tensor on large scales, i.e. ]J<1, one 
can actually give a measure of the amplitude of the scalar metric perturbation from the observations of 
the Sachs-Wolfe plateau in the angular power spectrum [20], 
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n = 0.97 ±0.03. 



(237) 



These measurements can be used to normalize the primordial spectrum and determine the parameters 
of the model of inflation [81 ]. In the near future these parameters will be determined with much better 
accuracy, as described in Section 4.4.5. 



4.43 The acoustic peaks 

The Sachs-Wolfe plateau is a distinctive feature of Fig. 24. These observations confirm the existence of a 
primordial spectrum of scalar (density) perturbations on all scales, otherwise the power spectrum would 
have started from zero at I = 2. However, we see that the spectrum starts to rise around / = 20 towards 
the first acoustic peak, where the SW approximation breaks down and the above formulae are no longer 
valid. 

As mentioned above, the first peak in the photon distribution corresponds to overdensities that 
have undergone half an oscillation, that is, a compression, and appear at a scale associated with the size 
of the horizon at last scattering, about 1° projected in the sky today. Since photons scatter off baryons, 
they will also feel the acoustic wave and create a peak in the correlation function. The height of the peak 
is proportional to the amount of baryons: the larger the baryon content of the universe, the higher the 
peak. The position of the peak in the power spectrum depends on the geometrical size of the particle 
horizon at last scattering. Since photons travel along geodesies, the projected size of the causal horizon 
at decoupling depends on whether the universe is flat, open or closed. In a flat universe the geodesies 
are straight lines and, by looking at the angular scale of the first acoustic peak, we would be measuring 
the actual size of the horizon at last scattering. In an open universe, the geodesies are inward-curved 
trajectories, and therefore the projected size on the sky appears smaller. In this case, the first acoustic 
peak should occur at higher multipoles or smaller angular scales. On the other hand, for a closed universe, 
the first peak occurs at smaller multipoles or larger angular scales. The dependence of the position of 

— 1/2 

the first acoustic peak on the spatial curvature can be approximately given by Z pca k ^ 220 Q , where 
Qo = + = 1 — ? ast observations from the balloon experiment BOOMERANG [87|, 
suggested clearly a few years ago that the first peak was between I = 180 and 250 at 95% cl., with an 
amplitude 5T = 80 ± 10 //K, and therefore the universe was most probably flat. However, with the high 
precision WMAP data we can now pinpoint the spatial curvature to a few percent, 

ft = 1-02 ± 0.02 (95% cl.) (238) 

That is, the universe is spatially flat (i.e. Euclidean), within 2% uncertainty, which is much better than 
we could ever do before. 




Multipole I k [h/Mpc] 



Fig. 29: The dependence of CMB anisotropies and LSS power spectrum on the sum of the mass of all neutrino species. The 
blue(red) data corresponds to WMAP(Boomerang, etc.) and SDSS(2dFGRS), for the CMB and LSS respectively. 

With Boomerang, CBI, VSA, and specially with WMAP, we have evidence of at least three distinct 
acoustic peaks. In the near furture, even before Planck, we may be able to distinguish anothes two. These 
peaks should occur at harmonics of the first one, but are typically much lower because of Silk damping. 
Since the amplitude and position of the primary and secondary peaks are directly determined by the 



sound speed (and, hence, the equation of state) and by the geometry and expansion of the universe, 
they can be used as a powerful test of the density of baryons and dark matter, and other cosmological 
parameters. 

By looking at these patterns in the anisotropics of the microwave background, cosmologists can de- 
termine not only the cosmological parameters, but also the primordial spectrum of density perturbations 
produced during inflation. It turns out that the observed temperature anisotropics are compatible with 
a scale-invariant spectrum, see Eq. d237l) . as predicted by inflation. This is remarkable, and gives very 
strong support to the idea that inflation may indeed be responsible for both the CMB anisotropics and the 
large-scale structure of the universe. Different models of inflation have different specific predictions for 
the fine details associated with the spectrum generated during inflation. It is these minute differences that 
will allow cosmologists to differentiate between alternative models of inflation and discard those that do 
not agree with observations. However, most importantly, perhaps, the pattern of anisotropics predicted 
by inflation is completely different from those predicted by alternative models of structure formation, 
like cosmic defects: strings, vortices, textures, etc. These are complicated networks of energy density 
concentrations left over from an early universe phase transition, analogous to the defects formed in the 
laboratory in certain kinds of liquid crystals when they go through a phase transition. The cosmolog- 
ical defects have spectral properties very different from those generated by inflation. That is why it is 
so important to launch more sensitive instruments, and with better angular resolution, to determine the 
properties of the CMB anisotropics. 

4.44 The new microwave anisotropy satellites, WMAP and Planck 

The large amount of information encoded in the anisotropics of the microwave background is the reason 
why both NASA and the European Space Agency have decided to launch two independent satellites to 
measure the CMB temperature and polarization anisotropies to unprecendented accuracy. The Wilkinson 
Microwave Anisotropy Probe [88] was launched by NASA at the end of 2000, and has fulfilled most of 
our expectation, while Planck [89 1 is expected to be lanched by ESA in 2007. There are at the moment 
other large proposals like CMB Pol [95 1, ACT l96l . etc. which will see the light in the next few years, 
see Ref. CHI . 

As we have emphasized before, the fact that these anisotropies have such a small amplitude allow 
for an accurate calculation of the predicted anisotropies in linear perturbation theory. A particular cosmo- 
logical model is characterized by a dozen or so parameters: the rate of expansion, the spatial curvature, 
the baryon content, the cold dark matter and neutrino contribution, the cosmological constant (vacuum 
energy), the reionization parameter (optical depth to the last scattering surface), and various primordial 
spectrum parameters like the amplitude and tilt of the adiabatic and isocurvature spectra, the amount of 
gravitational waves, non-Gaussian effects, etc. All these parameters can now be fed into very fast CMB 
codes called CMB FAST [93 1 and CAMB [94|, that compute the predicted temperature and polarization 
anisotropies to better than 1% accuracy, and thus can be used to compare with observations. 

These two satellites will improve both the sensitivity, down to //K, and the resolution, down to arc 
minutes, with respect to the previous COBE satellite, thanks to large numbers of microwave horns of var- 
ious sizes, positioned at specific angles, and also thanks to recent advances in detector technology, with 
high electron mobility transistor amplifiers (HEMTs) for frequencies below 100 GHz and bolometers for 
higher frequencies. The primary advantage of HEMTs is their ease of use and speed, with a typical sen- 
sitivity of 0.5 mKs 1 / 2 , while the advantage of bolometers is their tremendous sensitivity, better than 0.1 
mKs 1 / 2 , see Ref. [97 1. This will allow cosmologists to extract information from around 3000 multipoles! 
Since most of the cosmological parameters have specific signatures in the height and position of the first 
few acoustic peaks, the higher the resolution, the more peaks one is expected to see, and thus the better 
the accuracy with which one will be able to measure those parameters, see Table 2. 

Although the satellite probes were designed for the accurate measurement of the CMB tempera- 
ture anisotropies, there are other experiments, like balloon-borne and ground interferometers [90 1. Prob- 



ably the most important objective of the future satellites (beyond WMAP) will be the measurement of 
the CMB polarization anisotropics, discovered by DASI in November 2002 [98], and confirmed a few 
months later by WMAP with greater accuracy l20l . see Fig. 24. These anisotropics were predicted by 
models of structure formation and indeed found at the level of microKelvin sensitivities, where the new 
satellites were aiming at. The complementary information contained in the polarization anisotropics 
already provides much more stringent constraints on the cosmological parameters than from the temper- 
ature anisotropics alone. However, in the future, Planck and CMB pol will have much better sensitivities. 
In particular, the curl-curl component of the polarization power spectra is nowadays the only means we 
have to determine the tensor (gravitational wave) contribution to the metric perturbations responsible for 
temperature anisotropics, see Fig. |30| If such a component is found, one could constraint very precisely 
the model of inflation from its spectral properties, specially the tilt [91 1. 
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Fig. 30: Theoretical predictions for the four non-zero CMB temperature-polarization spectra as a function of multipole moment, 
together with the expectations from Planck. From Ref. 1 92 1. 



4.5 From metric perturbations to large scale structure 

If inflation is responsible for the metric perturbations that gave rise to the temperature anisotropies ob- 
served in the microwave background, then the primordial spectrum of density inhomogeneities induced 
by the same metric perturbations should also be responsible for the present large scale structure [99|. 
This simple connection allows for more stringent tests on the inflationary paradigm for the generation 
of metric perturbations, since it relates the large scales (of order the present horizon) with the smallest 
scales (on galaxy scales). This provides a very large lever arm for the determination of primordial spectra 
parameters like the tilt, the nature of the perturbations, whether adiabatic or isocurvature, the geometry 
of the universe, as well as its matter and energy content, whether CDM, HDM or mixed CHDM. 

4.51 The galaxy power spectrum 

As metric perturbations enter the causal horizon during the radiation or matter era, they create density 
fluctuations via gravitational attraction of the potential wells. The density contrast 5 can be deduced from 
the Einstein equations in linear perturbation theory, see Eq. dl66t . 

Sk ^ = (±-\*U i = (±-\*^l lk , (239) 
k p \aHj 3 * \aHJ 5 + 3w fe ' V ' 



where we have assumed K = 0, and used Eq. (I198I) . From this expression one can compute the power 
spectrum, at horizon crossing, of matter density perturbations induced by inflation, see Eq. (I210t . 



P(k) = (\5 k \ 2 ) = A (^)\ (240) 

with n given by the scalar tilt (I212L n = 1 + 2rj — 6e. This spectrum reduces to a Harrison-Zel'dovich 
spectrum (llOOt in the slow-roll approximation: ij, t « 1. 

Since perturbations evolve after entering the horizon, the power spectrum will not remain con- 
stant. For scales entering the horizon well after matter domination (k^ 1 S> k~ q x ~ 81 Mpc), the metric 
perturbation has not changed significantly, so that 72& (final) = 1Z k (initial). Then Eq. ( 12391 ) determines 
the final density contrast in terms of the initial one. On smaller scales, there is a linear transfer function 
T(k), which may be defined as [76 1 

72* (final) = T{k) ^(initial) . (241) 

To calculate the transfer function one has to specify the initial condition with the relative abundance 
of photons, neutrinos, baryons and cold dark matter long before horizon crossing. The most natural 
condition is that the abundances of all particle species are uniform on comoving hypersurfaces (with 
constant total energy density). This is called the adiabatic condition, because entropy is conserved inde- 
pendently for each particle species X, i.e. 5px = Px^t, given a perturbation in time from a comoving 
hypersurface, so 

_JP^_ = (242) 

PX +PX PY + PY 

where we have used the energy conservation equation for each species, px = —3H(px + px), valid to 
first order in perturbations. It follows that each species of radiation has a common density contrast S r , 
and each species of matter has also a common density contrast S m , with the relation 5 m = jS r . 

Given the adiabatic condition, the transfer function is determined by the physical processes oc- 
curing between horizon entry and matter domination. If the radiation behaves like a perfect fluid, its 
density perturbation oscillates during this era, with decreasing amplitude. The matter density contrast 
living in this background does not grow appreciably before matter domination because it has negligible 
self-gravity. The transfer function is therefore given roughly by, see Eq. (11031) . 

[ 1 , k <§C k e q 

T(k) = < 2 (243) 



The perfect fluid description of the radiation is far from being correct after horizon entry, because 
roughly half of the radiation consists of neutrinos whose perturbation rapidly disappears through free 
streeming. The photons are also not a perfect fluid because they diffuse significantly, for scales below 
the Silk scale, kg 1 ~ 1 Mpc. One might then consider the opposite assumption, that the radiation has 
zero perturbation after horizon entry. Then the matter density perturbation evolves according to 

5 k + 2H5 k + (c 2 s kl h -47TGp)5 k = 0, (244) 

which corresponds to the equation of a damped harmonic oscillator. The zero-frequency oscillator defines 
the Jeans wavenumber, kj = AirGp/cj. For k <C kj, 5 k grows exponentially on the dynamical 
timescale, Td yn = Imcj -1 = (AnGp)^ 1 ^ 2 = r grav , which is the time scale for gravitational collapse. 
One can also define the Jeans length, 




(245) 



which separates gravitationally stable from unstable modes. If we define the pressure response timescale 
as the size of the perturbation over the sound speed, r pre s ~ A/c s , then, if r pres > r grav , gravitational 
collapse of a perturbation can occur before pressure forces can response to restore hydrostatic equilibrium 
(this occurs for A > Aj). On the other hand, if r pres < r grav , radiation pressure prevents gravitational 
collapse and there are damped acoustic oscillations (for A < A j). 

We will consider now the behaviour of modes within the horizon during the transition from the 
radiation (c 2 = 1/3) to the matter era (c 2 = 0). The growing mode solution increases only by a factor of 
2 between horizon entry and the epoch when matter starts to dominate, i.e. y = 1. The transfer function is 
therefore again roughly given by Eq. (1243b . Since the radiation consists roughly half of neutrinos, which 
free streem, and half of photons, which either form a perfect fluid or just diffuse, neither the perfect 
fluid nor the free-streeming approximation looks very sensible. A more precise calculation is needed, 
including: neutrino free streeming around the epoch of horizon entry; the diffusion of photons around 
the same time, for scales below Silk scale; the diffusion of baryons along with the photons, and the 
establishment after matter domination of a common matter density contrast, as the baryons fall into the 
potential wells of cold dark matter. All these effects apply separately, to first order in the perturbations, to 
each Fourier component, so that a linear transfer function is produced. There are several parametrizations 
in the literature, but the one which is more widely used is that of Ref. [ 1 00 1 , 
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We see that the behaviour estimated in Eq. (1243 1 is roughly correct, although the break at k = k cq is not 
at all sharp, see Fig. |^ The transfer function, which encodes the soltion to linear equations, ceases to 
be valid when the density contrast becomes of order 1. After that, the highly nonlinear phenomenon of 
gravitational collapse takes place, see Fig.OTI 
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Fig. 31: The CDM power spectrum P(k) as a function of wavenumber k, in logarithmic scale, normalized to the local abun- 
dance of galaxy clusters, for an Einstein-de Sitter universe with h = 0.5. The solid (dashed) curve shows the linear (non-linear) 
power spectrum. While the linear power spectrum falls off like fc~ 3 , the non-linear power-spectrum illustrates the increased 
power on small scales due to non-linear effects, at the expense of the large-scale structures. From Ref. 1411 . 



4.52 The new redshift catalogs, 2dF and Sloan Digital Sky Survey 

Our view of the large-scale distribution of luminous objects in the universe has changed dramatically 
during the last 25 years: from the simple pre- 1975 picture of a distribution of field and cluster galaxies, 
to the discovery of the first single superstructures and voids, to the most recent results showing an almost 
regular web-like network of interconnected clusters, filaments and walls, separating huge nearly empty 
volumes. The increased efficiency of redshift surveys, made possible by the development of spectro- 
graphs and - specially in the last decade - by an enormous increase in multiplexing gain (i.e. the ability 
to collect spectra of several galaxies at once, thanks to fibre-optic spectrographs), has allowed us not 
only to do cartography of the nearby universe, but also to statistically characterize some of its properties, 
see Ref. [ 101 ]. At the same time, advances in theoretical modeling of the development of structure, with 
large high-resolution gravitational simulations coupled to a deeper yet limited understanding of how to 
form galaxies within the dark matter halos, have provided a more realistic connection of the models to the 
observable quantities 1 1 02 1 . Despite the large uncertainties that still exist, this has transformed the study 
of cosmology and large-scale structure into a truly quantitative science, where theory and observations 
can progress side by side. 

I will concentrate on two of the new catalogs, which are taking data at the moment and which 
have changed the field, the 2-degree-Field (2dF) Catalog and the Sloan Digital Sky Survey (SDSS). The 
advantages of multi-object fibre spectroscopy have been pushed to the extreme with the construction of 
the 2dF spectrograph for the prime focus of the Anglo- Australian Telescope E2l . This instrument is 
able to accommodate 400 automatically positioned fibres over a 2 degree in diameter field. This implies 
a density of fibres on the sky of approximately 130 deg -2 , and an optimal match to the galaxy counts 
for a magnitude bj ~ 19.5, similar to that of previous surveys like the ESP, with the difference that 
with such an area yield, the same number of redshifts as in the ESP survey can be collected in about 10 
exposures, or slightly more than one night of telescope time with typical 1 hour exposures. This is the 
basis of the 2dF galaxy redshift survey. Its goal is to measure redshifts for more than 250,000 galaxies 
with bj < 19.5. In addition, a faint redshift survey of 10,000 galaxies brighter than R = 21 will be done 
over selected fields within the two main strips of the South and North Galactic Caps. The survey has 
now finished, with a quarter of a million redshifts. The final result can be seen in Ref. \42\. 

The most ambitious and comprehensive galaxy survey currently in progress is without any doubt 
the Sloan Digital Sky Survey 1431 . The aim of the project is, first of all, to observe photometrically the 
whole Northern Galactic Cap, 30° away from the galactic plane (about 10 4 deg 2 ) in five bands, at limiting 
magnitudes from 20.8 to 23.3. The expectation is to detect around 50 million galaxies and around 10 8 
star-like sources. This has already led to the discovery of several high-redshift (z > 4) quasars, including 
the highest-redshift quasar known, at z = 5.0, see Ref. 1431 . Using two fibre spectrographs carrying 320 
fibres each, the spectroscopic part of the survey will then collect spectra from about 10 6 galaxies with 
r' < 18 and 10 5 AGNs with r' < 19. It will also select a sample of about 10 5 red luminous galaxies with 
r' < 19.5, which will be observed spectroscopically, providing a nearly volume-limited sample of early- 
type galaxies with a median redshift of z ~ 0.5, that will be extremely valuable to study the evolution of 
clustering. The data that is coming from these catalogs is so outstanding that already cosmologists are 
using them for the determination of the cosmological parameters of the standard model of cosmology. 
The main outcome of these catalogs is the linear power spectrum of matter fluctuations that give rise to 
galaxies, and clusters of galaxies. It covers from the large scales of order Gigaparsecs, the realm of the 
unvirialised superclusters, to the small scales of hundreds of kiloparsecs, where the Lyman-a systems 
can help reconstruct the linear power spectrum, since they are less sensitive to the nonlinear growth of 
perturbations. 

As often happens in particle physics, not always are observations from a single experiment suffi- 
cient to isolate and determine the precise value of the parameters of the standard model. We mentioned 
in the previous Section that some of the cosmological parameters created similar effects in the tem- 
perature anisotropies of the microwave background. We say that these parameters are degenerate with 



respect to the observations. However, often one finds combinations of various experiments/observations 
which break the degeneracy, for example by depending on a different combination of parameters. This is 
precisely the case with the cosmological parameters, as measured by a combination of large-scale struc- 
ture observations, microwave background anisotropics, Supernovae la observations and Hubble Space 
Telescope measurements. It is expected that in the near future we will be able to determine the param- 
eters of the standard cosmological model with great precision from a combination of several different 
experiments. 

5. CONCLUSION 

In the last five years we have seen a true revolution in the quality and quantity of cosmological data 
that has allowed cosmologists to determine most of the cosmological parameters with a few percent 
accuracy and thus fix a Standard Model of Cosmology. The art of measuring the cosmos has developed 
so rapidly and efficiently that one may be temped of renaming this science as Cosmonomy, leaving the 
word Cosmology for the theories of the Early Universe. In summary, we now know that the stuff we are 
made of — baryons — constitutes just about 4% of all the matter/energy in the Universe, while 25% is 
dark matter — perhaps a new particle species related to theories beyond the Standard Model of Particle 
Physics — , and the largest fraction, 70%, some form of diffuse tension also known as dark energy — 
perhaps a cosmological constant. The rest, about 1%, could be in the form of massive neutrinos. 

Nowadays, a host of observations — from CMB anisotropics and large scale structure to the age 
and the acceleration of the universe — all converge towards these values, see Fig. 25. Fortunately, we 
will have, within this decade, new satellite experiments like Planck, CMBpol, SNAP as well as deep 
galaxy catalogs from Earth, to complement and precisely pin down the values of the Standard Model 
cosmological parameters below the percent level, see Table 1. 

All these observations would not make much sense without the encompassing picture of the infla- 
tionary paradigm that determines the homogeneous and isotropic background on top of which it imprints 
an approximately scale invariant gaussian spectrum of adiabatic fluctuations. At present all observations 
are consistent with the predictions of inflation and hopefully in the near future we may have information, 
from the polarization anisotropics of the microwave background, about the scale of inflation, and thus 
about the physics responsible for the early universe dynamics. 
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